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SUMMARY 


A Continuously  Operating  Fast  Flow  Electrically  Excited 
(COFFEE)  laser  excitation  concept  was  studied  to  determine  the  feasibility 
of  scaling  it  to  high-power  levels.  The  concept  uses  a multipin-to-plane 
electrode  configuration  with  transverse  flow  to  produce  a high-pressure, 
sel*~sustained  glow  discharge.  Simultaneous  approaches  involving  both 
experimental  measurements  and  theoretical  analyses  were  employed  to 
develop  discharge  scaling  criteria  and  a laser  kinetics  model.  The 
results  demonstrate  that  the  COFFE^  laser  is  scalable  to  high-power  levels 
and,  in  addition,  exhibits  many  useful  and  unique  features. 

Experimental  measurements  were  made  of  the  glow  behavior  using 

an  electrode  module  (10  x 10  x 10  cm).  Parametric  studies  determined  the 

influences  of  discharge  current,  gas  velocity,  pressure,  mixture,  electrode 

gap,  and  cathode  pin  spacing  on  input  power  density,  specific  input  energy, 

discharge  stability,  and  scale  size.  Measurements  of  the  small-signal 

gain  distributions  were  also  made  for  a number  of  parameter  settings. 

In  addition  the  kinetic  gas  temperature  was  inferred  from  measurements 

of  the  small-signal  gain  taken  at  various  P and  R-branch  laser  wavelengths. 

The  results  included  observations  of  specific  input  energies  up  to  140  kJ/lb 

with  input  power  densities  up  to  50  W/cm3  of  discharge  volume.  Maximum 

iiput  power  was  obtained  at  ^300  Torr  for  a number  of  different  operating 

conditions.  With  an  8 cm  gap  a typical  small-signal  gain  of  0.25%/cm 

3 

was  obtained  at  30  W/cm  input  power  and  150  m/sec  gas  velocity.  The 
corresponding  gas  mixture  was  1C02: 7N2:20He:0.2H2  at  280  Torr.  These 
results  were  compared  with  the  theoretical  predictions  of  a glow  discharge 
model  aid  a laser  kinetics  model. 

The  glow  discharge  was  modeled  by  calculating  the  steady  state 
balance  between  electron  generation  and  loss  mechanisms  under  self-sustained 
operating  conditions.  This  model  predicted  the  operating  E/N  in  the 


*C/ 


upstream  portion  ot  the  discharge  to  within  10%  if  one  assumed  that  the 
electron  loss  rate  is  attachment,  dominated. 

The  laser  kinetics  were  analyzed  using  a computer  code  developed 
for  a six-temperature  Boltzmann  equilibrium  model  of  CO 7 :N  ; :He :H  0 gas 
mixtures.  Ihe  model  is  quite  general  and  can  be  applied  to  most  CO^ 
flowing  gas  laser  systems  of  present  interest.  Gas  temperature  and  E/N 
dependencies  were  included,  along  with  all  important  sources  of  gas  heating. 
Computer  predictions  of  the  small-signal  gain  distribution  along  the  gas 
flow  direction  agree  well  with  the  measured  values.  Laser  cavity  optics 
were  also  included  in  the  model,  and  parametric  calculations  were  performed 
for  various  output  couplings  and  operating  pressures.  An  output  power  of 
45  kW  at  an  efficiency  of  15%  was  predicted  from  an  8 x 11  x 140  cm  cavity 
operating  at  280  Torr. 

The  main  conclusion  of  these  studies  is  that  multikilowatt 
COFFEE  laser  systems  can  easily  be  ouilt  and  a solid  technology  base  now 
exists  for  designing  these  systems. 


1.  INTRODUCTION 


The  approach  in  this  study  program  is  based  on  the  use  of  a 
Continuously  Operating  Fast  Flow  Electrically  Excited  (COFFEE)  laser 
concept.  The  COFFEE  laser  uses  a new  high-pressure  (0.1  to  1 atm)  diiect 
current  self-sus'.ained  glow  discharge  excitation  scheme  which  is  well 
suited  for  a continuous  wave  C02  laser.  The  continuous  glow  is  produced 
with  a mult ipin-to-p lane  electrode  geometry  and  is  stabilized  by  a gas 
flow  transverse  to  the  discharge  and  by  a ballast  resistor  in  series 
with  each  cathode  pin.  The  excitation  concept  is  illustrated  in  Fig.  1:1. 
A high-power  CC>2  laser  system  using  this  excitation  scheme  would  be 
relatively  simple,  rugged,  compact,  and  inexpensive. 

The  objective  of  the  program  described  herein  has  been  the 
development  of  the  high-pressure  glow  discharge  scaling  criteria  for 
the  mult ipin-to-p lane  electrode  system.  Particular  emphasis  has  been 
placed  on  understanding  the  physics  of  the  discharge  as  they  relate 
to  discharge  scalability  and  laser  excitation  efficiency.  Parallel, 
complementary  approaches  involving  both  experimental  studies  and 
theoretical  modeling  were  pursued.  Specific  tasks  included  glow  discharge 
experiments  using  an  electrode  module  (10  x 10  x 10  cm)  to  measure  the 
important  glow  discharge  parameters.  This  module  was  constructed  during 
Phase  I of  this  study  program  as  described  earlier.1  The  results  for 
Phase  II,  which  was  a continuation  of  the  studies,  are  reported  in 
Section  2 of  this  report.  Small-signal  gain  measurements  are  also 
described.  Section  3 contains  the  theoretical  modeling  developed  to 
gain  further  understanding  of  the  discharge  characteristics  and  laser 
kinetics.  Preparatory  work  for  an  envisioned  Phase  III  to  build  a 25  kW 
proof  of  scalability  laser  is  described  in  Section  4.  Finally,  Section  5 
is  devoted  to  suggestions  for  further  work. 
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Fig.  1:1  Artist's  conception  of  COFFEE  laser 


2 . GLOW  DISCHARGE  EXPERIMENTS 


2.1  DESCRIPTION  OF  GLOW  DISCHARGE  SCHEME 


The  COFFEE  laser  concept  uses  a continuous,  self -sustained 
glow  discharge  which  is  produced  with  a multipin-to-plane  electrode 
geometry  as  shown  in  Fig.  1:1.  The  glow  is  stabilized  at  high  pressures 
by  a gas  flow  transverse  to  the  discharge  and  by  a ballast  resistor  in 
series  with  each  cathode  pin.  Additional  row  ballast  resistors  are  used 
to  tailor  the  input  power  distribution  as  a function  of  downstream  distance. 

As  part  of  the  COFFEE  laser  study  program,  an  electrode  module 
(10  x 10  x 10  cm)  was  tested  to  determine  the  glow  discharge  scalability 
as  a function  of  the  important  parameters.  One  purpose  of  the  module 
experiments  was  to  demonstrate  that  indeed  it  is  feasible  to  build  a 
large-scale  laser  device.  In  addition,  parametric  studies  were  made  to 
determine  the  influence  of  the  main  parameters,  such  as  discharge  current, 
gas  velocity,  pressure,  mixture,  electrode  gap,  and  cathode  pin  spacing. 
Small-signal  gain  measurements  were  made  as  a function  of  these  variables 
to  determine  the  uniformity  and  efficiency  of  the  excitation.  These 
experiments  have  also  provided  important  information  for  the  theoretical 
models  of  the  glow  discharge  and  laser  kinetics. 

The  discharge  module  formed  a two-dimensional  flow  nozzle  (sub- 
sonic) which  consisted  of  a planar  anode  on  one  side  and  a sheet  of  high- 
temperature  dielectric  on  the  other  side.  The  two  electrodes  assemblies  are  shown 
in  Figs.  2:1  and  2:2.  The  cathode  consisted  of  a multipin  array  (110  pins) 


protruding  approximately  0.6  cm  through  the  insulating  sheet  with  a pin 
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density  of  1 pin/cm  . Looking  along  the  10  cm  optical  axis,  transverse 
to  the  flow,  one  saw  11  rows  of  10  pins.  Each  pin  in  a given  row  was 
connected  through  a separate  50  LT  ballast  resistor  to  a common  row  bus  bar. 
As  the  laser  gas  mixture  passes  from  row  to  row,  the  discharge  impedance 
decreases  due  to  gas  heating.  This  was  compensated  in  the  present  case  by 
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adding  a separate  ballast  resistor  to  each  row  in  order  to  maintain  a 

nearly  constant  input  power  density  A diverging  flow  channel  could  also 
have  been  used. 

The  electrode  gap  was  varied  from  2 to  10  cm.  The  distance 
between  the  first  and  last  rows  of  cathode  pins  in  the  gas  flow  direction 
was  10  cm.  Alternate  rows  w.  : staggered  so  that  the  nearest  nin  in  the 

gas  flow  direction  was  2 cm  away,  minimizing  the  discharge  interactions 
at  the  cathode.  Tungsten  pins  with  a 1 mm  diameter  and  a sharp  conical 
tip  were  used.  The  copper  anode  was  15  cm  wide  and  approximately  30  cm 
long.  An  internally  metalized  glass  tube  was  fitted  around  the  perimeter 
of  the  anode  (not  shown)  to  act  as  guard  ring  for  minimizing  the  field 
enhancement  of  the  anode  edges  at  large  electrode  spacings.  The  internal 
metal  surface  of  this  tube  was  connected  electrically  to  the  downstream 
end  of  the  anode.  The  electrodes  were  fitted  into  a special  vacuum-tight 
fiberglass  test  section  housing  as  shown  in  Fig.  2:3.  This  test  section 
was  part  of  a closed-loop  wind  tunnel  capable  of  operating  at  any  pressure 
up  to  1 atm  with  variable  flow  rates  up  to  170  m3/min. 

2,2  ELECTRODE  MODULE  (10  x 10  x 10  cm)  PARAMETERIZATION 

The  high-pressure  direct  current  self-sustained  glow  discharge 
has  been  successfully  scaled  to  a 10  x 10  cm  square  aperture.  A photograph 
of  the  glow  discharge  is  shown  in  Fig.  2:4.  Visually  the  glow  appeared 
fairly  uniform,  although  no  direct  correlation  was  found  between  the 
measured  small-signal  gain  and  the  luminosity.  The  glow  vas  less  visible 
under  the  first  two  rows  because  the  current  to  these  rows  was  less  than 
to  the  other  rows.  The  upstream  rows  acted  to  "condition"  or  preionize  the 
gas  for  the  other  rows.  In  addition,  the  glow  discharge  appeared  to  be 
deflected  downstream  with  respect  to  the  cathode  by  approximately  15-30 
degrees,  depending  on  the  flow  velocity.  (See  Section  2.3  on  gain 
measurements).  For  spacings  up  to  10  cm  there  was  no  problem  with  the  glow 
blowing  out  or  becoming  positionally  unstable  for  gas  velocities  up  to 
150  m/sec.  The  discharge  luminosity  was  observed  to  decrease  again  under 
the  most  downstream  rows  where  the  current  density  was  the  highest.  This 
may  be  due  to  a slightly  lower  operating  E/N,  ratio  of  electric  field 
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Flow  150  m/sec 


High-pressure  DC  glow  discharge  viewed  along 
optical  axis  of  10  x 10  x 10  cm  nodule.  140  Torr 
lCO„:7N.:20He. 
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strength  to  gas  density,  in  this  region  as  described  in  detail  in 
Section  3.2.  Since  the  mechanisms  which  produce  luminosity  are  different 
from  those  which  produce  laser  gain,  one  should  not  expect  an  exact 
correlation  between  them. 

Average  input  power  densities  up  to  50  W/cm  of  discharge  volume 
were  obtained  with  specific  input  energies  up  to  140  kJ/lb  (>150  J/i-atm) . 

^ ^ iarge  gaps,  however,  somewhat  lower  values  were  obtained.  For  example, 
with  an  8 cm  discharge  gjw  we  obtained  30  W/cm  uniformly,  except  for 
the  first  two  cathode  rows,  at  50  kJ/lb  using  a ICO., : 7N9  : 201Ie  gas  mixture 
at  286  Torr  and  150  m/sec  flow  velocity.  This  result  is  shown  in 
Fig.  2.5  for  row  ballast  resistor  configuration  //l.  The  input  power 
distribution  can  easily  be  varied  over  a wide  range  as  demonstrated  by 
the  results  for  row  ballast  configuration  // 2.  Both  examples  in  Fig.  2:5 
corresponded  to  operating  conditions  just  below  the  glow-to-arc  threshold, 
and  the  total  input  power  turned  out  to  be  nearly  the  same  in  both  cases. 
Therefore,  no  attempt  was  made  to  optimize  the  row  ballast  configuration 
for  each  test  condition  during  the  discharge  parameterization  studies. 
Rather,  a near  optimum  configuration  for  uniform  input  power  density  at 
maximum  total  input  power  was  chosen  and  then  held  constant  during  each 
specific  test. 

The  influence  of  gas  velocity  on  the  maximum  average  input 
power  density  is  illustrated  in  Fig.  2:6  for  a typical  set  of  operating 
conditions.  In  general,  the  higher  velocities  tend  to  inhibit  the 
propagation  of  discharge  streamers  across  the  gap  and,  thereby,  prevent 
the  onset  of  arcing.  The  reason  for  the  nearly  linear  dependence  is 
poorly  understood  at  this  time.  As  described  previously,^  the  glow-to- 
arc  transition  appears  to  be  due  to  local  thermal  instabilities  in  the 
discharge  rather  than  bulk  thermal  effects.  Numerous  investigators^  ^ 
have  studied  the  phenomena  and  have  concluded  that  the  instability 
threshold  appears  to  be  related  to  many  factors,  such  as  gas  residence 
time,  input  power  density,  gas  turbulence  level,  and  local  field 
perturbations.  However,  the  present  theories  are  mainly  qualitative 
and  can  only  predict  general  trends.  Much  additional  work  is  needed  in 
this  area.  Since  the  highest  performance  was  desired,  most  of  the  other 
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Fig.  2:5  Variation  of  input  power  density  with  downstream 
distance  in  the  discharge  for  two  row  ballast 
resistor  configurations. 
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parameterizations  were  done  with  ^150  ra/sec  flow  velocity  which  was  the 
maximum  achievable  with  the  pressure  head  capability  of  our  blower. 

Figure  2:7  shows  the  maximum  total  input  power  to  the  discharge 
module  that  was  obtained  without  arcing  as  a function  of  discharge  gap 
for  two  operating  pressures.  The  solid  curves  are  straight  line 
approximations  of  the  experimental  data.  Higher  input  powers  were 
obtained  with  the  larger  gaps.  The  scalability  is  further  indicated  by 
replotting  this  data  in  terms  of  the  average  input  power  density  as  shown 
in  Fig.  2:8.  The  input  power  density  was  calculated  on  the  basis  of 
dividing  the  total  input  power  by  the  discharge  current  area  of  .1  cm  per 
pin  times  110  pins  (from  the  cathode  pin  density)  and  the  discharge  gap 
(2  to  10  cm).  The  solid  curve  is  a cross-plot  of  the  straight  line  from 
Fig.  2:7.  If  the  maximum  input  power  increased  linearly  with  gap,  the 
maximum  input  power  density  would  be  independent  of  the  discharge  gap, 
which  is  a desirable  scaling  relationship.  At  larger  gaps  (5-10  cm)  we 
see  that  the  maximum  input  power  density  was  fairly  insensitive  to  gap. 

At  small  gaps  (^2  cm),  however,  it  was  possible  to  achieve  almost  double 
the  input  power  density.  The  explanation  for  this  behavior  is  rooted  in 
the  discharge  instability  phenomena.  The  discharges  between  real  small 
gaps  are  probably  more  stable  because  the  component  of  the  electric  field 
which  presumably  exists  in  the  gas  flow  direction  is  considerably  smaller 
in  this  situation.  Nevertheless,  the  results  of  Fig.  2:7  indicate  that 
gaps  greater  than  10  cm  would  be  useful  for  obtaining  even  greater  power 
levels . 

The  effect  of  gas  pressure  is  shown  most  clearly  in  Fig.  2:9. 
The  average  input  power  densities  were  observed  to  increase  linearly 
with  discharge  current  for  any  given  operating  pressure  as  expected, 
because  the  V I characteristic  for  this  glow  discharge  is  nearly  flat 
(slightly  rising).  For  any  given  current  setting  the  input  power  density 
is  approximately  proportional  to  the  total  gas  pressure  if  the  gas 
mixture  is  held  cons  cant.  This  occurs  because  the  operating  E/N  for  the 
discharge  is  to  first  order  governed  only  by  the  gas  mixture.  The 
dashed  curve  in  Fig.  2:9  represents  the  glow-to-arc  limit  on  input  power 
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The  maximum  power  input  to  the  glow  discharge  without  sparking 
was  found  to  be  proportional  to  the  cathode  pin  array  area.  For  example, 
with  an  8 cm  gap  at  286  Torr  of  1C02 : 7^:  20He:  0. 6H2  mixture  and  flow 
velocity  of  150  m/sec,  it  was  possible  to  obtain  12.6  kW  input  using  only 
the  5 upstream  rows  of  cathode  pins.  With  all  11  rows  of  cathode  pins 
conoacj-*d,  r,.t  maximum  ^rc-ftct  input  power  was  ai.9  kW  (/uJOW/cc,  see 
Fig.  2:5).  A "quicky”  experiment  demonstrated  that  it  is  possible  to 
extend  the  discharge  downstream  to  more  than  20  cm  by  increasing  the 
cathode  array  dimension  in  the  flow  direction.  The  only  restriction, 
other  than  arcing,  is  that  one  must  not  overheat  the  gas  to  the  point 
that  lasing  action  stops  as  described  in  Section  3. 

The  pin  density  was  varied  from  0.25  to  1.00  pins/cm2  with 
negligible  affect  on  the  glow  stability  maximum  input  power.  This  variation 
was  accomplished  by  disconnecting  intermediate  rows  of  pins.  The  discharge 
was  not  as  uniform  in  the  vicinity  of  the  cathode  in  these  cases.  The  pin 
density  was  also  constrained  by  the  maximum  current  that  can  be  drawn  from  a 
single  pin.  For  example,  at  p = 280  Torr  with  a 1:7:20  gas  mixture 
flowing  at  150  ra/sec,  100  mA  was  about  the  maximum  single  pin  current 
that  could  be  obtained  with  a 4 cm  gap  and  upstream  discharges.  When 
only  a single  row  of  pins  having  the  same  1 cm  transverse  spacing  was 
inserted  into  this  same  gas  flow  field,  35  mA/pin  was  the  maximum  attainable 
current  withouu  arcing.  This  further  indicates  that  the  upstream  discharges 
tend  to  favorably  "condition"  the  gas  for  the  downstream  discharges. 

If  one  increases  the  pin  density,  there  is  still  an  upper 
limit  to  the  maximum  current  density,  based  on  the  area  of  the  discharge 
column,  which  can  be  obtained  before  arcing  occurs.  Adding  more  cathode 
Pins  just  lowers  the  current  per  pin  required  to  achieve  this  threshold 
current  density  and  makes  the  cathode  more  complex  to  build. 


p 


2-3  SMALL-SIGNAL  GAIN  MEASUREMENTS 


^ ^ • 3-  Description  of  Gain  Measurement  Technique 

The  experimental  setup  is  schematically  shown  in  Fig.  2:10. 

A tunable  C02  probe  laser  was  operated  in  the  TEM  fundamental  mode  to 
obtain  minimum  beam  divergence.  An  adjustable  iris  attenuated  the  beam 
to  intensity  levels  which  avoided  saturation  of  the  detector.  Germanium 
beam  splitters,  placed  along  the  path  of  the  laser  beam,  directed  samples 
of  the  beam  to  a power  meter,  a spectrum  analyzer  and  a UV  illuminated 
phosphor  screen.  The  laser  beam  was  then  chopped  with  a polished  stainless 
steel  blade  so  that  the  beam  reflected  off  the  chopper  passed  through 
a variable  attenuator  and  fell  on  a pyroelectric  detector.  The  transmitted 
beam  was  directed  through  the  laser  discharge  by  means  of  a set  of  beam 
steering  mirrors.  Any  region  of  the  discharge  along  a path  perpendicular 
to  the  gas  flow  direction  and  parallel  to  the  electrode  surface  could  be 
probed  by  moving  these  mirrors  horizontally  and  vertically.  A front 
surface  mirror  was  located  on  the  other  side  of  the  discharge.  Its 
alignment  was  adjusted  with  micrometer  screws  so  that  the  probing  beam 
was  reflected  almost  back  on  itself.  Thus,  the  returned  beam  also  fell 
on  the  same  pyroelectric  detector  as  shown  in  Fig.  2:10.  The  output  of 
the  pyroelectiic  detector  was  fed  into  a lock-in  amplifier.  The  reference 
signal  for  the  lock-in  amplifier  was  generated  by  chopping  a light 
source  with  the  same  chopper  mentioned  above.  It  was  detected  by  a 

photodiode.  The  output  of  the  lock-in  amplifier  was  displayed  on  a 
chart  recorder. 

Tne  system  was  initialized  by  first  observing  output  V from 
the  lockr-in  amplifier  when  the  beam  reflected  off  the  chopper  was 
temporarily  fully  attenuated.  Care  was  taken  to  assure  that  the 
pyroelectric  detector  was  operated  in  its  linear  range.  Then  the 
variable  attenuator  was  readjusted  until  the  lock-in  amplifier  output 
was  zero  with  the  main  C02  discharge  turned  off.  Finally,  the  discharge 
was  turned  on,  and  the  amplification  of  the  probe  beam  produced  an  output 
V2  from  the  lock-in  amplifier.  Thus,  the  total  measured  gain  was  given 
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by  V2^V1'  Since  the  total  active  path  length  was  20  cm,  the  percent 


gain  per  centimeter  was 


V, 


So  20  ^V1')  X 100 


(2:1) 


In  principle,  very  accurate  measurements  can  be  obtained  with  this  method 
because  extremely  small  values  of  V can  be  measured  using  this  nulling 
technique  and  the  high  sensitivity  scales  of  the  lock-in  amplifier.  The 
main  limitation  is  the  stability  of  the  reference  signal  V^.  In  the 
present  case  the  probe  laser  had  an  inherent  drift  which  limited  the  useful 
sensitivity  enhancement  to  a factor  of  ten.  This  permitted  the  measurement 
of  gain  as  low  as  0.01%/cm. 


2.3.2  Experimental  Results 

The  integrated  small— signal  gain  over  a 20  cm  path  perpendicular 
to  both  the  gas  flow  direction  and  the  electric  field  was  investigated  for 
a cathode-to-anode  spacing  of  8 cm.  The  variation  of  the  gain  in  the 
gas  flow  direction  (horizontal) , midway,  between  the  cathode  and  the  anode 
was  measured.  A vertical  scan  of  the  gain  was  also  made  between  the 
electrodes  at  the  position  where  maximum  gain  was  observed  in  the  horizontal 
scan.  The  dependence  of  the  gain  at  the  ?^2  transition  upon  gas  pressure, 
discharge  current,  and  hydrogen  partial  pressure  were  studied  at  the  mid- 
plane. The  gain  variation  with  J-values  of  both  the  P-branch  and  the 
R-branch  were  also  studied  for  a fixed  gas  composition  and  pressure. 

The  small-signal  laser  gain  for  the  ?22  trar-ition  as  a function  of  down- 
stream distance  from  the  first  row  of  cathode  pins  is  shown  in  Figs.  2:11 
and  2:12  for  total  discharge  currents  of  1.0  A and  1.5  A,  respectively. 

The  gas  mixture  was  ICO,,:  7N,, : 20He:  0. 2H,,.  The  total  gas  pressure  varied 
about  20%  during  the  course  of  the  measurements  because  some  gas  escaped 
out  of  the  blower  shaft  seal.  Hence,  the  gain  was  normalized  to  280  Torr 
using  the  pressure  dependence  shown  in  Fig.  2:14.  The  gas  flow  velocity 
was  150  m/sec,  and  the  specific  electrical  energy  input  to  the  discharge 
at  1.5  A was  50  kJ/lb.  The  row  ballast  resistors  were  adjusted  to  give 
nearly  constant  power  input  per  unit  volume  (see  Fig.  3:10).  At  1.5  A 
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the  discharge  was  operating  just  below  the  sparking  limit.  The 
variation  of  power  input  per  unit  volume  in  the  flow  direction  is 
further  described  in  Sections  2.2  and  3.2.4.  In  Figs.  2:11  and  2:12, 
the  dots  are  tht.  experimental  measurements  and  the  solid  curves  are  the 
theoretical  predictions  based  on  computer  calculations  with  various 
assumed  values  of  water  vapor  content  as  explained  in  Section  3.2.4. 

The  theoretical  curves  for  the  assumed  value  of  0.25%  water  content 
agrees  with  both  the  maximum  gain  and  the  shape  of  the  gain  curve 
within  the  scatter  of  the  experimental  data. 

In  order  to  match  the  position  where  the  experimental  and 
theoretical  curves  peaked,  the  theoretical  curve  was  displaced  downstream 
by  1.5  cm.  This  shift  is  qualitatively  consistent  with  a displacement  of 
the  current  distribution  arising  from  the  gas  flow.  The  drift  velocity 
of  positive  ions  in  the  discharge  is  given  by  W = u+  • E cm/sec  where  y 
is  the  mobility  of  positive  ions  in  the  mixture  and  E the  electric  field 


.in  volts/cm.  The  positive  ion  mobility  in  the  C0„:N„:He  mixture  used  is 
2 ^ 
about  7.1  cm  /volt-sec  and  the  E/p  under  the  last  row  of  cathode  pins 


was  6.1  volts/cm-Torr.  This  yields  a positive  ion  drift  velocity  of 


3.3  x 10  cm/sec.  The  gas  flow  velocity  here  was  1.6  x 10  4 cm/sec.  The 


current  path  appears  to  be  deflected  downstream  with  respect  to  the  cathode 
pins,  which  anchor  the  discharge,  by  an  angle  which  is  proportional  to  the 
ratio  of  gas  to  ion  drift  velocities.  Thus,  the  expected  displacement  of 
current  distribution  at  the  mid-plane  (4  cm  below  the  cathode)  was 
approximately  2 cm,  which  is  in  agreement  with  the  displacement  required 
to  match  theory  with  experiment. 

Figures  2:11  and  2:12  show  that  the  small-signal  gain  increased 
(by  integration)  as  the  laser  gas  passed  through  the  discharge.  Beyond 
the  discharge  region,  the  gain  fell  off  exponentially  with  a decay  time 
of  about  260  ysec.  This  relatively  long  decay  time  is  due  to  the 
nitrogen  acting  as  an  energy  reservoir.  The  shape  of  the  gain  profile 
in  the  discharge  is  governed  by  both  the  pumping  rate  and  the  deactivation 
rate.  The  details  of  these  processes  are  discussed  in  Section  3.  However, 
it  should  be  pointed  out  that  the  shape  of  the  gain  profile  in  the  down- 
stream direction  will  be  completely  different  in  the  presence  of  stimulated 
emission  during  lasing. 
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A vertical  a can  of  the  smalt-signal  gam  tor  the  P10  line  at 
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1.5  cm  downstream  from  the  last  row  of  cathode  pins  is  shown  in  Fig. 

2:13.  The  gain  expressed  in  percent  per  centimeter  is  plotted  against 
distance  measured  from  the  cathode.  Several  measurements  were  taken  at 
1.0  and  1.5  A at  each  vertical  position  and  were  normalized  to  a total 
pressure  of  280  Torr  using  a ICO.,:  7N, : 20He:  0.  2H.  fe-as  mixture. 
of  these  measurements  for  1.0  A are  indicated  by  dots  and  those  for  1.5  A 
are  indicated  by  triangles.  The  error  bars  correspond  to  the  standard  deviation 
of  the  data  taken  at  each  position.  The  results  indicate  that  within 
the  accuracy  of  the  measurements  the  gain  was  quite  uniform  across  the 
8 cm  interelectrode  spacing,  except  perhaps  for  a slight  fall-off  near 
the  cathode.  The  pressure  dependence  of  the  small-signal  gain  at  the 
P22  line,  for  a discharge  current  of  1.0  A and  for  a fixed  gas  composition 
of  1C02: 7N2: 20He:0. 2H2  in  the  pressure  range  of  70  to  280  Torr  is  shown 
in  Fig.  2:14.  The  measurements  were  made  midway  between  the  cathode  and 
the  anode  and  1 cm  downstream  from  the  last  row  of  cathode  pins.  The 
gas  flow  velocity  at  the  entrance  to  the  discharge  region  was  150  m/sec. 


The  log-log  plot  of  the  variation  of  the  gain  with  pressure  P shows  that 

,-.74 


the  gain  varies  asymptotically  as  P^,/‘*  for  large  pressures.  A similar 
pressure  dependence  was  also  observed  for  a 1.5  A discharge. 

Introduction  of  a small  amount  of  hydrogen  into  the  laser  gas 
mixture  improved  the  long  term  discharge  stability  conditions.  However, 
it  also  resulted  in  reduced  gain.  Fig.  2:15  shows  the  variation  of 
measured  small-signal  gain  with  total  discharge  current  for  various  amounts 
of  hydrogen  added  to  a fixed  1C02 : 7N2 : 20He  composition  at  a total  pressure 
of  ^280  Torr.  The  measurements  were  made  at  the  midplane  1 cm  downstream 
from  the  last  row  of  cathode  pins.  The  solid  curves  are  straight  lines 
fitted  through  the  data  points.  The  gain  varies  linearly  with  current 
for  a fixed  gas  composition  as  expected.  The  1.5  A data  are  replotted  in 
Fig.  2:16  as  a function  of  hydrogen  concentration.  In  general,  the 
small-signal  gain  was  reduced  approximately  40%  with  the  addition  of  1% 

H2  to  the  1:7:20  gas  mixture  at  280  Torr.  The  gain  decreased  exponentially 
with  further  addition  of  H2 . The  dashed  line  represents  the  theoretical 
predictions  for  various  assumed  values  of  water  vapor  content.  Unfortunately, 
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no  measurements  were  made  with  H20  added  in  place  of  H^,  A linear 


extrapolation  of  the  data  back  to  zero  H2  concentration  gives  a 
gain  equivalent  to  that  predicted  by  a 0.15%  H20  content.  One  might 
reasonably  expect  a flow  system  background  water  vapor  concentration 
of  this  magnitude  since  commercial  grade  gases  were  used  and  the 
system  was  not  bakable  under  vacuum.  The  reduction  in  gain  caused  by 
the  addition  of  H2  was  1/3  times  that  predicted  for  H20.  This  large 
reduction  in  the  measured  gain  implied  that  in  the  presence  of  the 
discharge,  chemical  reactions  occurred  which  led  to  the  formatio^,  of 
additional  water  vapor.  A more  detailed  comparison  of  the  two  curves 
is  given  in  Section  3.2.4. 


2 . 4 GAS  TEMPERATURE  MEASUREMENTS 


The  kinetic  temperature  in  molecular  gas  laser  media  such  as 


CC>2  can  be  inferred  from  measurements  of  the  small-signal  gain  taken 


at  various  P and  R-branch  laser  wavelengths.  This  technique  has  been 

, 5 


applied  by  Avizonis,  Dean  aid  Grotbeck.  using  P-branch  gain  measurements 
performed  upon  CC>2  gas  dynamic  lasers.  More  recently,  the  gain  formulations 
have  been  modified  at  Westinghouse  and  applied  to  P and  R-branch  gain 
measurements  in  UV-initiated,  pulsed  CCyN^He  laser  discharges.6  These 
gain  techniques  have  the  advantage  that  the  gas  temperature  is  determined 
at  the  site  and  time  of  the  laser  population  inversion,  and  the  method 
does  not  perturb  the  gain  medium.  For  the  COFFEE  laser  discharges  studied 
here,  P and  R-branch  gain  measurements  were  performed  in  the  interelectrode 
raidplane  1.0  cm  downstream  from  the  eleventh  row  of  pi.s  for  total  currents 
of  1.0  and  1.5  A. 


The  small-signal  gain  relationship  at  line  center  can  be 
expressed  as^ 


.qQlA(J)  = k b ,n  , exp  [-J(J  ± 1)  6 /T  ] 
J F(v, J)  c v v 


" KcBv"V'  exp  + 1)6v"/Tr-1 


(2:2) 
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for  a vibrational-rotational  transition  from  the  (v',J')  level  to  the 

(v",J")  level.  The  lower  signs  are  selected  for  P-branch  gain  with  the 

rotational  quantum  number  J = J"  designated  by  the  lower-lying  level. 

The  R-branch  gain  is  obtained  by  selecting  the  upper  signs  with  J = J' 

designated  by  the  higher-lying  level.  The  terms  a(J),  A(J)  and  F(v,J) 

are  the  measured  small-signal  gain,  the  wavelength  corresponding  to  the 

particular  transition,  and  the  vibrational-rotational  interaction  factor, 

respectively.  The  quantities  and  n^  are  the  rotational  constant  and 

the  population  density,  respectively,  for  the  v^  vibrational  level.  The 

characteristic  rotational  temperature  6^  is  defined  in  terms  of  the 

rotational  constant  B as 

v 

Gv=hcBv/k,  (2:3) 

where  h,  c and  k are  Planck's  constant,  the  velocity  of  light  and 
Boltzmann's  constant,  respect.'. vely . 

The  equilibrium  rotational  temperature  T^  is  assumed  to  be 
identical  for  both  the  upper  and  the  lower  laser  levels,  and  T^  is 
further  assumed  to  be  equal  to  the  gas  kinetic  temperature  T.  These 
assumptions  are  well-founded  for  the  CO^  lasers  considered  here,  since 
rotational  thermalization  and  equilibrium  with  the  gas  kinetic  reservoir 
occur  in  times  on  the  order  of  1 nsec.  Thus,  a determination  of  the 
rotational  temperature  T^  from  gain  measurements  provides  an  excellent 
estimate  of  the  gas  temperature  based  upon  the  assumption  that  T % T . 

The  collisional  constant  applies  for  collisional  broadening 
of  the  vibrational-rotational  transition,  and  is  given  by 

3 

K = 1671  oc  S , „/3kT  v (2:4) 

c v'v"  r c ' ' 


k 

Note:  This  definition  of  J for  the  R-branch  simplifies  the  equations 

in  this  section;  however,  the  equations  in  Appendix  E of  Section  3.1 
are  similarly  simplified  by  taking  the  alternate  convention  of  setting 
J = J"  for  the  R-branch. 
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where  a is  the  molecular  symmetry  number  (0  = 2 for  C0„),  and  Sv,v„  is 
the  square  of  the  vibrational  transition  moment.  The  collision  frequency 


f°r  optical  broadening  can  be  expressed  as 


v 


= (8kT/M)1/2£n.u.  (1  + M/M.) 1/2 

j J J J 


(2:5) 


where  M is  the  mass  of  the  laser  particle  and  n.,  0.  and  M are  the 

J J j 


number  density,  optical  broadening  cross  section  and  mass,  respectively, 
of  the  gas  mixture  constituents.  The  sum  is  taken  over  all  gas  species, 


including  the  laser  species.  Under  the  assumption  that  T = T,  the 

collisional  constant  K is  seen  to  vary  as  T -2//2 

c J r 


Inspection  of  Eq.  (2:2)  reveals  that  if  the  small-signal  gain 
u is  measured  for  various  transitions  (v',J')  ->  (v",J"),  the  only 
undetermined  quantities  are  nyf  , nv„  and  Tf.  Since  ny,  and  nv„  are  not 
known  a griori,  further  approximations  to  Eq . (2:2)  are  required  in  order 


to  extract  Tf  from  the  gain  measurements.  This  can  be  accomplished  by 


removing  the  first  term  of  Eq.  (2:2)  as  a common  factor  and  approximating 
the  remaining  term  in  exp[+2J8/Tr]  by  unity  with  Bv,  = Bv„  = B and 

6v'  = 6v"  = 0-  For  typical  C02  values  of  J = 20  with  Tf  = 300°K,  this 
exponential  term  differs  from  unity  by  only  about  7%;  thus  the 
approximation  does  not  introduce  significant  error  into  the  gain 
formulation  for  typical  C02  laser  conditions.  The  resulting  gain 
equation  is  considerably  simplified,  however: 


q(J)  MJ) 


J F(v,J) 


KcB(nv*  “ nv'^  exPl-JU  ± 1)  6/T^J 


(2:6) 


Taking  natural  logarithms  of  both  sides  of  this  equation  vields  a 
convenient  linear  relationship: 


£n[a(J)  X(J)/J  F(v, J)  ] = -J(J  + 1)  |-  + £n [K  B(nv,  - n „)] 


(2:7) 
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The  left-hand  side  involves  the  measurable  quantities  o(J)  and 
J,  and  the  first  term  on  the  right-hand  side  involves  the  measurable 
quantity  J(J  + 1).  Thus  a plot  of  £n[ct(J)  A(J)/J  F(v,J)]  as  a function 
of  J(J  + 1)  will  yield  the  same  straight  line  for  both  P and  R branch 
data.  The  slope  of  this  line  is  -6/T^  and  the  ordinate  intercept  is 

*n[KcB(V  " V)]-  The  §as  temperature  T = T^  can  be  evaluated  from 
the  slope,  and  the  population  difference  (nyl  - d^,)  can  be  cal^Ia.ed 
from  the  ordinate  intercept  value  and  T 

This  method  of  determining  the  gas  temperature  has  been  programmed 
for  automatic  computer  treatment  of  raw  experimental  gain  data.  The 
computer  code  GAHED  calculates  the  appropriate  variables  for  Eq.  (2:7), 
and  performs  a linear  least  squares  fit  of  these  values  to  a straight 
line.  Both  P-branch  and  R-branch  gain  data  are  considered  separately, 
since  the  approximation  exp(+2Je/Tj  1 provides  a slightly  different 
linear  relationship  for  P and  R-b ranches.  The  approximation  causes  an 
overestimation  of  Tr  from  the  P-branch  data,  and  an  underestimation  of 
Tr  from  the  R-branch  data.  Thus  separate  fitting  procedures  applied  to 
P and  R-branch  gain  data  should  bracket  the  true  value  of  T„;  in  practice, 
the  arithmetic  average  of  these  values  is  taken  as  the  gas  kinetic 
temperature. 

The  numerical  constants  employed  in  these  computations  were 
derived  from  the  latest  measured  values  reported  in  the  scientific 
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literature.  Wavelengths  for  the  various  CO,,  P and  R-branch  transitions 
were  taken  from  the  measurements  of  Barker  and  Adel,8  and  the 


vibrational-rotational  interaction  factors  F(v,J)  were  deduced  from  the 

f nrmiil  - ■ i . O 


num  uie 

formulations  given  by  Aria,  Lacome  and  Rossetti.9  Rotational  constants 

- — 0 . r*m  D r\  onm  “"1  _ 10 


"00*1  ’ °-3866  “d  Bio»o  ■ °-3897  from  Hertberg10  were  used 

the  C0_  UDDP.r  1 nijor  1 i _ . 
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for  the  C02  upper  and  lower  laser  levels,  respectively.  The  molecular 

cinnmft  ♦-  vu  i _ o ^ 


symmetry  number  a = 2 for  CO,,  was  obtained  from  Herzberg.11 


The  line 


-O'  ■‘•lie  lllie 

strength  S for  the  <00°1)  - (10°0)  laser  transition  in  CO,  was  taken 
as  1.48"?  y in  — d c . , 2 


__29  2 z. 

as  1.483  x 10  erg  cm  from  the  recent  measurements  of  Murray,  Kruger 


22  WJ-  ) ivLuge 

and  Mitchner,  and  the  CO,  optical  broadening  cross  sections  for  CO  , 
N-  flnH  Ho  Morn  10-1  1 -I  / . , _ 2 
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and  He  were  taken  as  131,  114  and  48  °AZ,  respectively  from  this  SL 

Or* on  r>  a 


ref erenc 


32 


T 


it 


i . 


>■  r; 


Gain  measurements  were  performed  at  a location  1.0  cm  downstream 
from  the  eleventh  row  of  pins  in  the  COFFEE  electrode  array,  midway  between 
the  anode  and  cathode.  An  8 cm  electrode  gap  was  used.  The  laser  gas 
mixture  was  1C02: 7N2: 20He: 0. 2H2  with  a total  pressure  of  280  Torr.  A 


nxed  entering  gas  velocity  of  150  m/sec  was  used.  With  a total  discharge 
current  of  1.0  A,  small-signal  gain  measurements  were  performed  at  8 
P-branch  and  7 R-branch  transitions.  At  the  1.5  A discharge  current, 

8 P-branch  and  10  R-branch  gain  measurements  were  made.  The  particular 


transitions  probed  were  carefully  chosen  to  avoid  the  anomalous  CO 
'hot  bands"  reported  by  Leonard13  and  by  Singer.14  Specifically,  no 


gain  measurements  were  attempted  on  the  P(20),  P(28)  and  R(15)  transitions 
of  the  10.4  |Jm  C02  bands,  since  these  transition  wavelengths  apparently 
overlap  those  corresponding  to  other  C02  levels  which  can  be  highly 
populated  at  typical  operating  temperatures.  Such  "hot  band"  transitions 
can  lead  to  anomalously  high  observed  gain  at  these  wavelengths,  by  as 
much  as  10%  in  some  cases.  Since  such  anomalies  could  influence  the 
rotational  temperatures  deduced  from  these  measurements,  the  "hot  band" 
wavelengths  were  not  included  in  the  data  set. 

The  observed  gains  for  the  various  P and  R-branch  C02  transitions 
are  shown  in  Figs.  2:17  and  2:18.  These  values  have  been  normalized  to  a 
constant  background  pressure  of  280  Torr  by  employing  the  empirically- 
determined  relationship 


a « p 


-0.7 


(2:8> 


where  a is  the  observed  gain  and  p is  the  measured  total  gas  pressure. 
Since  changes  in  background  pressure  were  recorded  over  the  1 to  4 hour 
duration  of  a typical  experimental  run,  this  normalization  procedure  was 
deemed  necessary  to  correct  for  an  observable  systematic  error.  Further- 
more, each  data  point  is  the  average  of  between  3 and  10  separate 
observations.  Nonetheless, it  can  be  seen  that  considerable  scatter  exists 
m the  experimental  gain  data.  This  is  attributed  to  instabilities  in 
the  output  of  the  probe  laser  and  the  rather  small  signal-to-noise  ratio 
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obtained  in  the  detection  equipment.  The  computer  least-squares  fit  to 
the  data  is  indicated  approximately  by  the  solid  curve  drawn  through  the 
experimental  points. 

For  the  1.0  A case,  the  P and  R branch  temperatures  determined 
from  the  GARED  program  were  301°K  and  326°K,  respectively,  to  yield  an 
average  of  313°K.  For  the  1.5  A case,  these  values  were  289°K  and  298°K, 
respectively,  for  an  average  of  293°K.  These  results  are  anomalous  in 
that  the  higher  current  case  yields  a temperature  lower  than  the  lower 
current  case,  and  the  1.5  A data  yields  a temperature  below  the  'v300°K 
inlet  gas  temperature.  Consequently,  we  conclude  that  the  scatter  in  the 
experimental  data  is  sufficient  to  introduce  significant  errors  into  the 
temperature  determination.  Based  upon  error  calculations  performed 
previously  for  UV-initiated,  pulsed  CO^  laser  gain  measurements,  we 
estimate  that  the  uncertainty  in  the  derived  temperatures  is  ^20°K 
for  the  1.0  A data  and  ^40°K  for  the  1.5  A data.  Clearly  an  improved  gain 
measurement  technique  will  be  required  to  reduce  errors  in  the  gas 
temperature  determination. 

Despite  these  uncertainties,  the  gain  reduction  technique 
demonstrated  that  the  gas  temperature  increased  by  ^30°K  during 
traversal  through  the  COFFEE  excitation  region.  Furthermore,  this  result 
agreed  with  the  laser  kinetics  prediction  shown  in  Fig.  3:11.  With  this 
temperature  rise,  it  is  evident  that  the  CO2  laser  kinetics  were  not 
affected  seriously  by  the  deleterious  effects  of  elevated  gas:  temperatures 
which  occur  beyond  ^450°-500°K.  Thus  excitation  regions  extended  even 
further  in  the  downstream  dimension  should  result  in  increased  laser 
excitation  and  power  output. 
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3.  THEORETICAL  MODELING 


3.1  LASER  KINETICS  MODEL 


3.1.1  Introduction 


In  the  COFFEE  CO^  laser  kinetics  model  we  have  considered  five  general 
mechanisms:  (1)  the  e]ectrical  excitation  of  the  gaseous  medium,  (2)  the 

collisional  relaxation  of  energy  in  the  gas,  (3)  the  transport  of  thermal 
energy  by  gas  flow,  (4)  the  interaction  of  coherent  radiation  with  the 
C02  molecules,  and  (5)  the  feedback  of  the  coherent  radiation  by  end 
mirrors.  Considerable  work  has  been  done  in  these  areas  by  many  workers, 
particularly  in  conventional  low  pressure,  low  energy  CC>2  devices.  The 
object  of  developing  the  present  model  was  to  extend  previous  treatments 
to  include  continuous  operation,  the  temperature  dependence  of  the 
collisional  relaxation  processes,  the  electric  field  dependence  of  the 
electrical  excitation,  and  the  effects  of  water  vapor  and  gas  flow.  The 
basic  model  is  the  one  which  has  been  successful  in  simulating  pulsed 
CO ^ lasers.  However,  it  was  modified  to  include  the  gas  flow  properties 

and  to  model  the  continuous  operation  of  the  COFFEE  laser.  The  complete 
COFFEE  C02  laser  kinetics  model  is  described  in  the  following  sections. 

3 .1  . 2 General  Description  of  the  Model 

The  basic  premise  of  the  COFFEE  C02  laser  kinetics  model  is  that 
the  vibrational  modes  of  the  molecules  can  be  considered  as  completely 
separate  energy  reservoirs  which  interact  by  electronic  and  molecular 
collisions  and  by  photon  emission  or  absorption.  Each  mode  is  assumed  to 
be  in  local  Boltzmann  equilibrium  due  to  rapid  intramode  vibrational- 
vibrational  (W)  relaxation.  Thus,  the  energy  content  of  each  vibrational 
mode  is  characterized  by  its  own  Boltzmann  temperature  and  is  described 
by  its  own  rate  equation. 


The  fraction  of  the  input  electrical  power  which  excites  the 

various  vibrational  modes  and  electronic  levels  has  been  determined  bv 
16  * 
Lowke  et  al.  from  electron  energy  distribution  function  computations  in 

CC^  laser  mixtures.  Thesr  fractions  are  computed  as  a function  of  E/N, 

where  E is  the  electric  field  and  N is  the  total  number  density  of  the  gas 

In  general,  the  value  of  E/N  in  the  discharge  is  a function  of  position 

downstream  and  its  value  at  any  position  was  determined  from  the  value 

of  the  electric  field  calculated  from  the  measured  glow  voltage  and  from 

the  total  number  density  calculated  from  the  specific  heat  rate  equation, 

assuming  that  the  pressure  is  constant  throughout  the  laser  cavity.  The 

experimental  values  of  the  current  density  were  used  to  obtain  the  input 


electrical  power. 

From  the  theoretical  models  one  is  able  to  predict  the  proper 
temperature  dependence  of  the  vibrational-vibrational  (W)  and  vibrational- 
translational  (VT)  collisional  relaxation  rate  constants.  However,  existing 
theories  yield  quite  poor  values  for  the  magnitude  of  these  rates.  Consequer 
the  theoretical  models  were  used  to  provide  a three-parameter  expression 
for  the  temperature  dependence  of  these  rate  constants,  and  this  expression 
was  fitted  to  measured  values  at  three  different  temperatures.  Selaxation 
rates  for  which  no  experimental  measurements  are  available  were  estimated 
on  the  basis  of  physical  arguments.  Most  of  the  rates  involving  water 
must  be  estimated  in  this  fashion.  Water  was  included  in  this  model 
because  its  relaxation  rates  are  extremely  high  and  exhibit  anomlous 
temperature  behavior.  Thus,  even  very  small  amounts  of  water  vapor  will 
have  a significant  effect  on  the  results.  On  the  other  hand,  hydrogen 
has  not  been  included  because  its  relaxation  rates  are  at  least  eight 
times  slower  than  those  of  water  vapor  and  because  the  concentration  of 
hydrogen  was  never  allowed  to  become  very  large. 

As  the  gas  temperature  increases  in  the  gas  flow  direction,  the 
gas  flow  velocity  increases  and  the  number  density  decreases.  These 
changes  are  directly  proportional  to  the  temperature  or  its  inverse  and 
were  included  in  the  model.  Although  these  changes  are  important,  the 
collisional  rates  are  also  dependent  on  the  gas  temperature  and  in  many 
cases  exhibit  large  changes  for  small  temperature  variations.  Thus,  it 
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is  necessary  to  consider  the  nonstationary  equation  of  thermal  conductivity 
to  determine  the  gas  temperature  as  a function  of  the  downstream  position 
(all  gas  properties  are  assumed  to  be  constant  within  any  cross-section 
perpendicular  to  the  gas  flow  direction).  The  model  employs  the  specific 
heats  of  the  gas  components  at  constant  pressure  and  includes  the  three 
major  sources  of  gas  heating:  energy  defect  heating  from  VV  collisional 
transitions,  intermode  heating  from  VT  collisional  transitions,  and  gas 
heating  from  the  electrical  discharge.  The  loss  of  electrical  input  energy 
by  electronic  excitation  followed  by  photon  emission  was  also  included. 

The  total  radiation  intensity  is  the  sum  of  the  spontaneous 
radiation  intensity  not  collimated  by  the  cavity  mirrors  and  the  spontaneous 
plus  stimulated  radiation  intensity  which  is  collimated  by  the  cavity 
mirrors.  The  collimated  radiation  intensity  was  calculated  from  the  rate 
equation  which  defines  the  laser  gain.  Included  were  the  effects  of 
spontaneous  emission,  constant  input  radiation  intensity  (in  case  the 
laser  is  used  as  an  amplifier),  mirror  and  cavity  losses.  Doppler,  natural, 
and  pressure  broadening  of  the  linewidth  were  also  included.  The 
fractional  populations  of  the  rotational  levels  were  considered  in  computing 
the  stimulated  emission  transition  rates. 

As  just  described,  the  computation  proceeds  as  shown  in 
Fig.  3:1.  The  experimental  values  of  the  electric  field,  E,  and  current 
density,  j , as  a function  of  downstream  position,  x,  were  used  to  determine 
the  electron  excitation  rate  coefficients  and,  the  input  power  values. 

This  information  was  used  in  the  laser  Kinetics  model  to  calculate  the  gas 
temperature,  T,  and  many  other  laser  parameters,  taking  into  consideration 
the  discharge  and  laser  cavity  dimensions  as  well  as  the  optics  design, 
initial  gas  flow  velocity,  temperature,  and  pressure.  The  collisional 
rate  coefficients  determine  the  collisional  excitation  and  relaxation.  The 
environmental  changes  a molecule  experienced  as  it  flowed  through  the  dis- 
charge cavity  were  incorporated  in  the  laser  kinetics  model.  The  result  was 
the  determination  of  the  gas  flow  velocity,  gas  density,  gas  temperature, 
vibrational  temperatures,  gain,  output  power,  and  output  efficiency  as  a 
function  of  the  downstream  position.  The  basic  rate  equations  used  in 
this  model  are  described  in  the  next  section. 


Fig.  3:1  Basic  flow  diagram  of  the  COFFEE  CO  2 laser 
kinetics  model  computations. 


3*1*3  The  Rate  Equations 


From  Che  basic  premise  of  separate  energy  reservoirs  for  each 
vibrational  mode,  the  rate  equations  can  be  derived  by  introducing  the 
following  approximations: 

1.  Translational  and  rotational  motions  are  so  tightly  coupled 
that  their  energy  contents  are  in  Boltzmann  equilibrium  with  each  other 
at  the  kinetic  gas  temperature  T. 

2.  Only  the  bending  mode  of  HO  (WBM)  interacts  with  the 
vibrational  levels  of  C02  and  N2;  the  other  two  ^0  vibrational  modes, 
lying  considerably  higher  in  energy,  do  not  exchange  significant  amounts 
of  vibrational  energy  with  C02  and  N?. 

3.  The  WBM  can  be  represented  as  a perfect  simple  harmonic 
oscillator  (SHO)  with  an  infinite  number  of  equally  spaced  energy  levels,1^ 
since  only  the  lowest  levels  of  the  H20  bending  mode  are  appreciably 
occupied  at  the  temperatures  of  interest. 

4.  The  N2  vibrational  mode  (NVM) , the.  CC>2  asymmetric  stretch 
mode  (ASM) , the  CC>2  symmetric  stretch  mode  (SSM)  , and  the  doubly 
degenerate  CC>2  bending  mode  (BM)  can  be  represented  as  perfect  simple 
harmonic  oscillators  with  an  infinite  number  of  equally  spaced  energy 
levels,  as  their  lower-lying  energy  levels  exhibit  nearly  uniform 
spacings. 

5.  The  intramode  W relaxation  is  so  rapid  that  each 
vibrational  mode  is  in  local  Boltzmann  equilibrium  and  can  be  characterized 
by  its  own  Boltzmann  temperature. 

6.  The  anharmonicities  in  the  ASM  and  SSM  are  large  enough  to 
shift  the  wavelengths  of  transitions  between  the  Onl  and  InO  energy  levels 
(n  •>  1)  away  from  the  00°1  to  10°0  transition  wavelength,  and  out  of  the 
wavelength  region  for  which  there  is  appreciable  overlap  of  the  linewidths. 
This  wavelength  shift  and  the  small  occupancies  of  the  upper  energy 

levels  permit  one  to  restrict  the  stimulated  emission  to  the  00°1  - 10°0 
transition. 

7.  Ihe  resonator  cavity  wavelength  discrimination  is  sufficient 
to  suppress  stimulated  emission  at  9.6  urn. 
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8.  Spontaneous  emission  and  particle  diffusion  rates  are 
sufficiently  small  that  these  processes  can  be  neglected.  (However, 
spontaneous  emission  is  included  in  the  gain  equation  to  provide  the 
initial  radiation  which  is  amplified  by  the  laser  whenever  it  is  used 
as  an  oscillator.) 

9.  The  gas  pressure  is  constant  in  time,  and  molecular 
dissociation  can  be  neglected. 

The  model  is  considerably  simplified  by  averaging  over  all 
spatial  dependencies  transverse  to  the  direction  of  flow  and  by  restricting 
the  laser  emission  to  a single  vibrational-rotational  transition  within 
the  10.4  pm  band.  The  former  simplification  restricts  the  resonator 
configuration  to  be  plane-parallel.  In  addition  to  these  approximations, 
it  is  assumed  that  transitions  into  or  out  of  the  doubly  degenerate 
bending  mode  occur  only  when  the  energy  transfer  involves  either  one 
or  the  other  bending  mode. 

It  is  further  assumed  that  the  application  of  the  Landau- 
Teller  selection  rules  for  vibrational  level  transitions  is  valid  for 
collisional  interactions  between  modes,  and  that  the  dependence  of 
transition  probabilities  on  vibrational  quantum  numbers  in  multiple  quanta 
exchanges  may  be  calculated  from  the  one-dimensional  exponential  repulsion 
model  of  Schwartz,  Slawsky,  and  Herzfeld.  For  brevity,  the  latter  model 
is  called  the  SSH  theory. 

The  specific  interactions  considered  in  the  model  are  illustrated 
in  Fig.  3:2  for  the  usual  CC^ :He gas  composition.  Included  in 
the  model  are  three  VT,  six  W,  one  stimulated  emission,  and  five  electrical 
excitation  processes.  The  model  also  included  all  inverse  transitions 
determined  from  the  principle  of  detailed  balance.  These  interactions 
include  all  those  known  to  be  important.  The  inclusion  of  water  vapor 

introduces  interactions  which  have  not  been  included  in  previous 
models.  The  stated  relaxation  times  are  typical  for  a 1:7:20  mixture  at 
280  Torr  pressure. 
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Fig.  3:2  Vibrational  modes  of  molecules  important 
CO„  laser  mixtures 
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The  energy  content  of  each  vibrational  mode  is  described  by  a 

temperature:  T for  the  NVM,  T for  the  WBM,  T for  the  ASM,  T for  the 
n w a s 

SSM,  and  for  the  BM.  The  interactions  in  Fig.  3:2  generate  a flow  of 
energy  between  the  vibrational  modes.  The  application  of  the  aforementioned 
assumptions  to  these  interactions  leads  to  a rate  equation  for  the 

21 

Boltzmann  factor  associated  with  each  of  the  five  vibrational  modes. 

The  inclusion  of  a rate  equation  for  the  kinetic  gas  temperature  T and 
one  for  the  coherent  light  intensity  completes  the  model.  The  result 
is  a harmonic  oscillator,  six-temperature  Boltzmann  equilibrium  model 
described  by  the  following  seven  simultaneous  differential  equations: 
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where 
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In  the  last  expression  v is  the  gas  fJo»  velocity  and  is  the  doTOStream 
distance  from  the  laser  cavity  entrance.  These  rate  equations  agree  in 


essence  with  those  derived  by  Lyon,  22  by  Granek, 23  and  by  Harrach  and 


Einwohner.  The  individual  symbols  are  explained  in  Appendix  A, 
Nomenclature-, 


The  collisional  relaxation  terms  contain  a factor  R for  VV 

“g 


relaxation  of  the  a mode  due  to  collisions  with  molecules  of  species  g. 
Similarly  the  factor  Ra  describes  VT  relaxation  of  the  a mode.  In 
Eq.  (3:3)  the  first  and  third  factors  of  2 arise  from  the  double  degeneracy 
of  the  C02  bending  mode.  The  factor  of  3 arises  from  the  exchange  of  one 
ASM  equation  for  three  BM  quanta,  while  the  second  factor  of  2 arises  from 
the  exchange  of  one  SSM  quantum  for  two  BM  quanta.  All  other  energy 
exchanges  involve  one  quantum  transitions.  Similar  reasons  explain  the 
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multiplicative  factors  of  2 and  3 in  Eq.  (3:6).  Each  collisional 
relaxation  rate  constant  was  given  by  a three-parameter  temperature- 
dependent  expression  which  was  fitted  to  the  measured  values  given  in 
the  work  of  Taylor  and  Bitterman,  "~>  or  estimated  on  the  basis  of  physical 
arguments.  The  details  of  this  parameterization  and  temperature  plots 
of  the  24  rate  constants  used  in  the  laser  kinetics  model  are  given  in 
Appendix  B,  Collisional  Relaxation. 

The  percentage  of  the  input  electrical  power  which  excites  the 

various  vibrational  modes  and  electronic  levels  was  calculated  from  the 

1 f) 

Boltzmann  equation  as  a function  of  E/l.  by  Lowke  et  al.  By  separating 

the  c Jinbined  vibrational  levels  according  to  their  degeneracies,  the 

results  given  in  Fig.  3:3  were  obtained  for  a 1:7:20:0.07  mixture  of 

C°2:N2:He  Excitation  of  the  water  vibrational  mode  at  concentrations 

below  1%  is  too  low  to  be  considered.  The  curves  indicate  that  an  E/N 

16  2 

value  of  about  2 x 10  V-cm  is  near  optimum  for  the  electrical 
excitation  of  the  nitrogen.  The  electrical  pumping  rates  of  the  various 
vibrational  modes  are  denoLed  by  in  the  rate  equations  and  are  easily 

derived  from  Fig.  3:3.  The  details  of  the  derivation  are  given  in 
Appendix  C,  Electrical  Excitation. 

The  gas  temperature  was  determined  from  the  nonstationary 
equation  of  thermal  conductivity,  Eq.  (3:6).  The  bracketed  term  on  the 
left-hand  side  of  the  equation  contains  the  specific  heats  of  the  gas 
components  at  constant  pressure  weighted  by  their  respective  moirr 
densities.  lhe  first  three  lines  on  the  right-hard  side  contain  the 
energy  defect  heating  terms  from  VV  transitions,  while  the  fourth  line 
contains  the  VT  heating  terms.  The  first  term  of  the  last  line  arises 
from  gas  heating  by  the  discharge.  This  term  is  modified  by  a factor 
which  accounts  for  the  portion  of  electronic  excitation  which  directly 
heats  the  gas  rather  than  relaxing  through  photon  emission.  The  following 
set  of  terms  in  parentheses  subtracts  that  portion  of  the  electrical 
discharge  energy  which  goes  directly  into  vibrational  excitation.  A 
detailed  explanation  of  this  equation  is  given  in  Appendix  D,  Gas 
Heating . 
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The  equations  and  terms  described  thus  far  are  sufficient  to 
accurately  describe  the  small-signal  gain.  However,  if  there  is 
significant  stimulated  emission  present,  the  stimu.lated  emission  rate 
F in  Eqs.  (3:1)  and  (3:2)  and  the  gain  Eq . (3:7)  musL  also  be  considered. 


The  latter  equation  describes  the  radiation  intensity  I'  collimated  by 

y3 


the  resonator  mirrors  where  the  subscript  specifies  the  particular 
vibrational-rotational  transition  being  considered.  This  radiation  is  the 
result  of  stimulated  emission  and  is  the  final  output  of  the  laser.  The 
first  term  on  the  right-hand  side  arises  from  the  net  gain  of  the  medium 
where  is  the  total  gain  and  is  the  averaged  loss  per  unit  length. 

The  averaged  loss  includes  all  mirror  losses,  including  mirror  transmissivities, 


and  scattering  losses  within  the  gas  medium.  This  loss  term  is  frequently 

26 

called  the  cavity  gain  and  is  the  gain  at  which  a steady  state  laser  will 


operate.  The  second  term  on  the  right-hand  side  acgounts  for  the  spontaneous 

emission  which  initiates  the  laser  process,  and  the  last  term  represents 

irradiation  of  the  laser  medium  by  an  external  source  of  intensity  I , 

o 

in  the  event  that  the  gain  medium  is  employed  as  an  amplifier  driven  by 
an  external  laser  oscillator.  The  stimulated  emission  rate  is  fully 
explained  in  Appendix  E,  Stimulated  Emission,  and  the  calculation  of  the 
radiation  intensity  is  described  in  Appendix  F,  Radiation  Intensity. 

The  harmonic  oscillator,  six-temperature  Boltzmann  equilibrium 
laser  kinetics  model  is  described  fully  by  these  seven  simultaneous 
differential  equations.  This  system  of  equations  was  integrated  numerically 
in  constant  integration  steps  by  the  Euler  method^  to  obtain  the  behavior 
of  the  laser  quantities  as  a function  of  their  downstream  location.  A 
typical  computer  calculation  took  40  minutes  on  the  U-1106  computer  for  a 
layer  cavity  20  cm  wide.  The  laser  cavity  was  considered  to  have  a 
rectangular  cross-section  and  to  be  at  least  as  large  as  the  discharge 
region  (see  Fig.  3:4).  The  mirrors  were  assumed  to  be  plane-parallel 
rectangular  mirrors  along  the  sides  of  the  discharge  region  and  large 
enough  to  assure  that  the  entire  discharge  region  was  part  of  the  resonant 
cavity.  These  assumptions  simplify  the  model  but  are  not  fundamental 
limitations.  The  solution  of  these  equations  yielded  considerable 
information  about  the  operation  of  the  C07  COFFEE  laser,  including 
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the  small-signal  gain  and  laser  output  power.  Examples  of  the  information 
obtained  from  these  calculations  and  how  this  information  has  enhanced 
our  understanding  of  the  COFFEE  laser  operation  are  given  in  the  next 
section . 

3.1.4  Results 

A good  method  for  testing  a laser  kinetics  model  and  computer 

code  is  to  compare  the  theoretical  calculations  with  the  measured  values 

for  typical  causes.  However,  care  musL  be  taken  to  include  significant 

amounts  of  all  the  gas  components  and  to  include  all  significant 

collisional  relaxation  processes.  The  case  must  also  be  in  a regime 

where  the  experiment  can  be  monitored  accurately  and  reliably.  The 

conditions  chosen  were  a 280  Torr,  1:7:20  laser  mixture  of  CO^iN^iHe 

with  2 Torr  of  in  a laser  cavity  11  cm.  along  the  gas  flow  direction, 

10  cm  long  (in  the  direction  between  the  mirrors) , and  having  an  8 cm 

gap  between  the  electrodes  (i.e.,  a laser  cavity  volume,  V,  of  about  1 l) . 

The  small-signal  gain  was  measured  for  these  conditions  for  a total 

electrode  current  of  1.0A  and  1.3A  with  entering  flow  velocities  of 

149  and  148  m/sec,  respectively. 

The  steady-state  electric  field  and  current  density 

distributions  measured  for  total  currenLs  of  1.0A  and  1.5A  are  shown  in 

Figs.  3:5  and  3:6.  The  interelectrode  current  density  associated  with 

each  cathode  pin  was  determined  by  measuring  the  current  to  each  pin 

2 

and  assuming  that  the  currenL  spreads  uniformly  to  occupy  a 1 cm  area 

which  corresponds  to  half  of  the  distance  to  the  next  pin  in  each 

2 

direction  (pin  density  = 1 pin/cra  ).  The  current  density  and  electric 
field  are  arbitrarily  taken  to  be  zero  0.5  cm  before  and  after  the 
discharge  region  as  defined  by  the  cathode.  This  cutoff  procedure  may 
be  slightly  too  abrupt,  but  it  shouldn't  affect  the  general  results. 

The  pin  currents  were  experimentally  adjusted  to  increase  with  downstream 
distance,  as  the  electric  field  decreased,  in  an  attempt  to  keep  a constant 
input  power  density  into  the  gas. 
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The  current;  density  and  electric  field  distributions  shown  in 
Figs.  3:5  and  3:6  are  for  the  midplane  between  the  cathode  pins  and  the 
anode  plate  where  the  small-signal  gain  was  measured.  These  distribution 
plots  were  purposely  shifted  1.5  cm  downstream  from  the  x origin 
(x  - 0 at  0.5  cm  upstream  from  the  first  row  of  pins)  in  order  to  match 
the  location  of  the  calculated  peak  gain  with  that  of  the  measured  peak 
gain  (see  Figs.  3:7  and  3:8).  The  relative  locations  of  the  rows  of  pins 
are  shown  in  the  figures.  The  shift  at  the  midplane  in  the  distributions 
can  be  explained  by  a vector  combination  of  the  ion  drift  velocity  and 
the  gas  flow  velocity.  A simple  calculation  predicts  a 2 cm  displacement 
in  the  flow  direction  of  the  glow  discharge  at  the  midplane  which  agrees 
reasonably  well  with  the  1.5  cm  shift  required  to  match  the  calculated 
and  measured  gain  distributions.  (See  Section  2.3.) 

The  measured  current  density  and  electric  field  distributions 
allow  two  cases  to  be  calculated:  the  small-signal  case  when  mirrors 

are  not  present  and  the  large-signal  case  when  mirrors  are  present.  In 
the  former  case  detailed  experimental  measurements  of  small-signal  gain 
were  taken  at  1.0A  and  1.5A  for  comparison  purposes.  However,  in  the 
latter  case  experimental  measurements  with  laser  output  were  not  made 
because  the  discharge  module  was  only  10  cm  long.  Nevertheless,  it  was 
informative  to  predict  the  laser  performance  using  the  current  density 
and  electric  field  distributions  measured  at  1.5A.  For  these  lasing 
calculations  the  assumed  cavity  length  was  140  cm  and  the  cavity  width 
and  electrode  gap  were  kept  at  11  cm  and  8 cm,  respectively,  giving  a 
cavity  volume  of  about  12  •.  and  total  currents  of  14A  and  21A  corresponding 
to  the  1.0A  and  1.5A  currents  of  the  smaller  cavity.  The  mirrors  were 
assumed  to  be  plane-parallel  sheet  mirrors  which  formed  the  two  ends  of 
the  cavity  (sides  of  flow  channel),  with  one  mirror  having  a transmissivity 
of  0%  and  the  other  15%.  The  absorption  and  scattering  losses  were  taken 
to  be  1.5%  at  each  mirror  and  0.01%  cm  ^ within  the  gaseous  discharge. 
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The  ?iuall-sig:ial  gain  results  £or  the  P(22)  transition  art  shown 
in  Figs.  3:7  and  3:8.  The  calculated  small-signal  gains  are  obviously 
too  high  if  no  water  vapor  is  included  in  the  gas  mixture,  whereas  they 
atree  very  well  with  the  experimental  values  wlvn  tfcu-  gas  mixture  is 
assumed  to  contain  0.25%  water  vapor,  i.e.,  a 1:7:20:0.07  gas  mixture 
of  CO^iN^iHe.'H^O.  The  theoretical  and  experimental  values  agree 
particularly  well  in  the  region  of  peak  gain  where  the  experimental  values 
are  most  accurate.  We  will  therefore  use  0.25%  water  vapor  in  the 
theoretical  calculations. 

In  Fig.  3:7  the  bump  in  the  gain  curve  for  x between  3 and  4 
cm  is  a result  of  the  dip  in  the  current  density  distribution  shown  in 
Fig.  3:5.  This  slight  dip  for  the  1A  case  resulted  from  using  the  row 
ballast  resistance  values  which  were  optimized  for  1.5A  operation. 

Typical  collisional  decay  times  for  the  1.5A  test  are  given  in  Fig.  3:2 
and  were  considerably  shorter  than  those  exhibited  in  Figs.  3:7  and  3:8. 
The  apparent  difference  between  these  results  is  easily  explained  by  an 
examination  of  the  basic  rate  Eqs.  (3 : 1—3 : 5) . The  decay  times  in 
Fig.  3:2  are  the  reciprocals  of  the  collisional  transition  rates  Ra  and 
g whereas  the  decay  times  exhibited  in  Figs.  3:7  and  3:8  are  the  result 
of  the  entire  terms  in  the  rate  equations,  i.e.,  they  include  the  effects 
of  back  rates  and  reduced  transition  probabilities  due  to  elevated 
vibrational  temperatures  which  increase  the  populations  of  the  upper 
energy  levels.  However,  the  decay  times  given  in  Fig.  3:2  effectively 
show  that  water  vapor  decreases  the  small-signal  gain  by  depopulating 
the  excited  nitrogen  more  effectively,  relative  to  other  depopulation 
mechanisms,  than  it  does  the  CO^  bending  mode. 

The  large-signal  gain  results  in  Figs.  3:7  and  3:8  follow  the 
small-signal  gain  until  the  onset  of  lasing  &t  6.5  cm  for  the  14A  case 
and  at  4.5  cm  for  the  21A  case).  At  that  location  the  gain  falls  to  a 
cavity  gain  value  of  0.074%  cm  ^ as  determined  by  the  assumed  cavity 
losses.  Finally,  at  x = 12.0  cm  the  gas  leaves  the  discharge  region; 
lasing  ceases,  and  the  remaining  gain  decreases  exponentially  due  to 
collisional  deactivation.  The  large-signal  gain  calculation  thus  predicts 

54 


Without  Mirrors 
With  Mirrors 
Experimental 
Measurements  without 
Mirrors 


Position  of  Cathode  Pins 


1 i 


laser  output  extending  from  x = 6.5  cm  to  12.0  cm  for  the  14A  case  and 
from  x = 4.5  cm  to  12.0  cm  for  the  21A  case.  The  shorter  lasing  interval 
for  the  14A  case  is  expected  from  the  lower  electrical  pumping  rate;  that 
is,  it  takes  longer  and,  therefore,  more  distance  to  achieve  the  same 
amount  of  pumping  as  for  a current  of  21A.  However,  the  shutoff  of  the 
lasing  radiation  generally  occurs  at  the  same  location  because  the  gain 
is  at  the  cavity  threshold  value  (gain  = losses)  irrespective  of  the 
current  value.  Any  small  decrease  below  this  threshold  gain,  in  this 
instance  due  to  rapid  fall-off  of  discharge  current,  means  that  the  gas 
becomes  primarily  on  absorptive  medium  for  the  lasing  radiation,  thereby 
quickly  dissipating  the  laser  radiation.  These  predicted  performance 
features  appear  to  be  quite  reasonable  and  provide  a good  basis  for 
guiding  future  lasing  experiments. 

Although  the  theoretical  and  experimental  small-signal  gain 
values  agree  well,  a discrepancy  is  apparent  when  the  two  gas  mixes  are 
compared.  The  total  pressure  was  280  Torr  in  both  cases;  however,  the 


theoretical  model  assumed  that  the  gas  mixture  contained  0.7  Torr  Ho0 


and  no  whereas  in  the  experiments  2.0  Torr  and  no  Ho0  were  added 


to  the  initial  gas  mixture.  This  discrepancy  may  be  better  understood 


by  comparing  the  peak  gains  calculated  as  a function  of  Ho0  concentration 


with  the  peak  gains  measured  as  a function  of  concentration  as  shown 


in  Fig.  3:9.  (For  brevity  only  the  1.5A  or  21A  case,  unless  stated 
otherwise,  will  be  considered  in  the  remainder  of  this  section.)  The 
experimental  value  at  0%  H0  has  not  been  included  in  the  straight  line 
fit  to  the  experimental  data  because  it  was  measured  with  fresh  gas 
(the  gas  had  not  yet  been  recycled  many  times)  in  a cool  system  whereas 
the  other  measurements  were  made  with  old  gas  (the  gas  had  been  recycled 
many  times)  in  a warm  system.  The  gain  of  fresh  gas  was  always  higher 


than  that  of  old  gas.  The  extrapolated  peak  gain  at  0%  li0,  as  shown  by 


the  dotted  lines,  corresponds  to  0.15%  H 0 0 . This  suggests  that  the 


experimental  gas  mixtures  always  contained  a background  impurity  Ho0 


concentration  of  1.5%,  a value  which  is  not  unreasonable.  However,  the 
calculated  and  measured  small-signal  gain  values  agreed  best  when  0.25% 
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H20  was  included  in  the  theoretical  model.  The  need  for  assuming  a 


larger  H20  concentration  can  be  explained  either  by  assuming  that  the 


collisional  relaxation  rates  of  H20  and  H2  are  simply  scalable,  or  by 
assuming  that  a chemical  reaction  occurs  whereby  some  of  the  H2  is 
converted  to  H90.  If  the  background  H20  concentration  is  0.15%,  Fig.  3:9 


shows  that  3.2  times  as  much  H2  is  needed  to  produce  the  same  effect 


as  a given  amount  of  H20.  On  the  other  hand,  a comparison  of  the 
measured  relaxation  rates25  shows  that  H20  is  8.0  times  more  effective 


than  H2  in  relaxing  the  C02  asymmetric  stretch  mode  to  the  (X>2  bending 
mode,  which,  compared  to  competing  mechanisms,  is  the  most  important 
process  for  the  small-signal  gain.  Although  the  numbers  are  admittedly 
approximate,  they  seem  to  indicate  that  the  results  in  Figs.  3:7  and 
3:8  cannot  be  achieved  by  a simple  scaling  of  the  H2  and  H20  relaxation 


rates. 


The  other  possible  explanation  is  that  part  of  the  H2  was 


chemically  converted  to  H20.  By  assuming  a background  of  0.15%  H 0, 


this  requires  that  an  additional  0.1%  H20  was  produced  by  chemical 


conversion.  This  implies  that  14%  of  the  H2  concentration  in  the 


experimental  gas  underwent  chemical  conversion  to  H20,  and  that  the 


remaining  H2,  because  of  its  lower  collisional  relaxation  rates,  had 


no  significant  effect  on  the  laser  kinetics.  A 14%  conversion  factor 
is  reasonable  and  gives  strong  support  to  the  hypothesis  that  important 
chemical  reactions  occur  in  the  gas  flowing  through  the  discharge. 

To  be  efficient  ^ laser  must  have  good  pumping  properties. 

One  such  desirable  property  is  to  have  a uniform  electrical  power  density 
deposited  into  the  gas  (due  to  gain  spikes  and  gas  temperature  effects 
downstream,  a constant  input  power  density  may  not  be  the  optimum  design, 
but  it  is  still  a good  design).  The  measured  power  density  distribution 
with  no  mirrors  is  shown  in  Fig.  3:10  for  the  1.5A  case  in  which  the  row 
ballast  resistors  were  chosen  to  produce  a uniform  input  power  density. 
The  resulting  j-E  distribution  was  very  uniform  in  the  downstream 
direction  through  the  optical  cavity.  This  illustrates  one  of  the 
unique  features  of  the  COFFEE  laser  electrical  excitation  concept. 
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Another  important  pumping  property  is  to  operate  at  the  value 
of  E/N  for  which  maximum  power  is  deposited  in  the  proper  molecular 
excitations.  Reference  to  Fig.  3:3  shows  that  for  a 1:7:20  gas  mixture, 


— 1 6 2 

the  optimum  E/N  value  for  the  N2  excitation  is  2.0  x 10  V-cm  . By 


using  the  experimental  values  of  electric  field  (see  Fig.  3:6)  and 
calculating  the  value  of  N as  the  gas  flows  downstream  and  its  temperature 
increases,  the  E/N  values  can  be  determined  as  a function  of  x as  shown  in 
Fig.  3:10.  The  COFFEE  laser  clearly  operates  very  near  the  optimum  value 
with  75  to  83%  of  the  electrical  power  being  deposited  in  the  N2 
vibrational  mode.  This  operating  condition  is  better  than  that  of  non- 
flowing C02:N2:He:H20  self-sustained  electrical  discharges  where  the 


1 fi  2 

predicted  E/N  value  must  be  a constant  3.0  x 10  V cm  with  68%  of  the 


electrical  power  being  deposited  in  the  N2  vibrational  mode. 


The  combination  of  uniform  power  density  deposition  and  excellent 
E/N  operating  values  means  that  the  COFFEE  laser  operates  at  near  optimum 
pumping  efficiencies  throughout  the  optical  cavity.  Although  the  currents 
and  electric  fields  were  measured  during  the  small-signal  gain  experiments, 
the  addition  of  mirrors  should  have  a minor  effect,  and  the  basic 
conclusion  of  near  optimum  pumping  efficiencies  should  remain  valid. 

The  gas  temperature  is  important  because  the  collisional 
relaxation  processes  are  strongly  dependent  upon  it.  The  temperature 
distributions  with  and  without  mirrors  are  shown  in  Fig.  3:11  where  a 
typical  starting  temperature  of  300°K  is  assumed.  (The  gas  flow  velocity 
is  proportional  to  the  gas  temperature  and  therefore  Fig.  3:11  is  also  a 
plot  of  the  flow  velocity  when  300°K  is  replaced  by  148  m/sec  and  340°K 
is  replaced  by  170  m/sec) . The  presence  of  mirrors  markedly  increases  the 
gas  temperature  in  the  region  of  electrical  excitation  between  x = 1.5 
and  12.0  cm  and  decreases  the  temperature  farther  downstream.  These 
features  are  easily  understood  by  examinimg  the  interactions  diagrammed 
in  Fig.  3:2.  In  the  small-signal  case  there  are  no  significant  10.6  pm 


lasing  transitions,  and  the  relaxation  of  the  C02  asymmetric  mode  to  the 


00^  bending  mode  is  by  slow  W interactions.  However,  in  the  large  signal 
case  there  are  significant  lasing  transitions  to  the  C02  symmetric  mode 
which  is  closely  coupled  to  the  C02  bending  mode.  Thus  in  the  electrical 
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excitation  region,  gas  heating  by  VT  relaxation  of  the  C02  bending  mode  is 
low  in  the  small-signal  case  and  high  in  the  large-signal  case.  In  other 
words,  the  increased  gas  temperature  is  caused  by  the  increased  VT 
relaxation  of  the  C02  bending  mode.  Eventually,  however,  the  gas  temperature 
in  the  small-signal  case  must  exceed  that  in  the  large-signal  case  since 

the  laser  radiation  energy  in  the  latter  case  is  not  available  to  heat  the 
gas. 


The  three  sources  of  gas  heating  are  VT  and  VV  relaxation  and 
direct  heating  by  the  electrical  discharge.  The  first  two  sources  occur 
by  molecule-molecule  collisions,  whereas  the  latter  source  occurs  by 
electron-molecule  collisions.  The  contributions  to  the  gas  heating  from 
these  three  sources  are  shown  in  Figs.  3:12  and  3:13.  In  both  case  the 
VV  heating  is  small  but  understandably  larger  in  the  small-signal  case, 
where  W relaxation  is  the  only  interaction  between  the  asymmetric  and 
bending  modes  of  CO^ . The  electrical  discharge  contribution  is  a function 
of  only  the  current  and  voltage  distributions,  and  is,  therefore,  the 
same  for  both  cases.  The  VI’  relaxation  collisions  give  the  largest 
contribution,  which  is  the  same  for  the  two  cases  until  the  laser  output 
is  initiated  at  x = 4.5  cm.  As  explained  in  the  previous  paragraph,  the 
laser  radiation  effectively  increases  the  CC>2  bending  mode  VT  relaxation, 
thereby  heating  the  gas.  This  effect  is  clearly  seen  in  these  two  figures. 

The  vibrational  mode  Boltzmann  temperatures  for  the  small  and 
large-signal  cases  are  plotted  in  Figs.  3:14  and  3:15.  As  expected,  the 
C02  asymmetric  mode  and  N2  vibrational  mode  are  in  close  equilibration  as 
are  the  H20  bending  mode  and  the  symmetric  stretch  and  bending  modes. 

In  fact,  the  latter  three  modes  are  in  such  close  equilibration  that 
their  temperatures  cannot  be  separated  in  the  figures.  In  Fig.  3:15  the 
absence  of  any  significant  vibrational  temperature  increase  in  the  region 
of  lasing  radiation  from  x = 4.5  to  12.0  cm  indicates  that  the  pumping 
power  deposited  into  the  N'2  vibrational  mode  and  the  C02  asymmetric  stretch 
mode  is  being  effectively  removed  by  the  laser  cavity  flux  before  it  can 
increase  the  Boltzmann  temperatures.  This  is  a clearly  desirable 
characteristic  of  a good  laser. 
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Flg.  3:12  Contributions  to  gas  heating  when  mirrors 
are  not  present. 


Curve  678697-A 


V = 121 
P = 280  Torr 
T = 300°K 
v = 148  m/sec 


Although  the  Boltzmann  temperatures  for  the  N2  vibrational  mode 
and  the  CO ^ asymmetric  stretch  mode  are  the  important  parameters  in  the 
laser  kinetics  model  and  are  related  to  the  amount  of  energy  sLored  in  the 
gas,  a more  directly  meaningful  quantity  is  Lhe  total  energy  density  stored 
in  these  two  modes  which  can  be  theoretically  extracted  at  10.6  pm,  assuming 
one  photon  for  each  vibrational  quantum.  This  quantity  is  plotted  in  Fig. 
3.16.  The  small-signal  case  shows  that  the  gas  is  still  being  effectively 
pumped  at  12  cm  although  not  as  efficiently  as  at  4 cm  where  the  slope  is 
higher.  The  large-signal  case  shows,  as  in  the  previous  figure,  that  in  the 
region  of  lasing  radiation , the  pumping  power  deposited  in  the  vibrational 
mode  and  in  the  CO^  asymmetric  stretch  mode  is  being  effectively  removed 
by  the  laser  cavity  flux  before  it  can  increase  the  stored  energy  density. 
From  this  figure  it  can  be  determined  that  14.1  kW  of  theoretically 
extractable  radiation  power  remains  in  the  gas  beyond  the  region  of 
electrical  excitation,  although  all  of  this  energy  can  never  be  extracted 
in  a real  laser. 

The  integrated  output  efficiency  and  power  are  shown  in  Fig. 

3:17.  The  final  output  power  of  45.4  kW  is  at  an  efficiency  of  15.2%  but 
the  efficiency  and  output  power  are  both  increasing  at  the  end  of  the 
eJectrical  excitation  region.  This  indicates  that  the  COFFEE  laser's 
output  power  and  efficiency  can  be  increased  by  simply  extending  the 
excitation  region  beyond  the  present  11  cm.  This  is  also  seen  in  Fig. 

3:11  which  shows  that  the  present  COFFEE  laser  is  not  temperature  limited 
since  the  335  K is  well  below  the  450°K  at  which  the  laser  performance 
starts  to  degrade  due  to  collisional  deactivation.  Furthermore,  the  14.1 
kW  of  theoretically  extractable  output  power  which  remains  in  the  gas 
past  the  excitation  region  (see  Fig.  3:16)  is  virtually  independent  rf 
the  excitation  region  size  and  is  difficult  to  recover  experimentally. 

Thus,  an  excitation  region  extended  in  the  flow  direction  lowers  the 
relative  power  loss  in  the  gas  leaving  the  excitation  region, thereby 
increasing  the  efficiency. 
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Curve  680387-A 


J 45 


^ / 12: 7: 20: 0.07 
/ V = 12  J0 
/ P = 280  Torr 
T = 300° K 
v = 148  m/sec 
t = 15% 
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The  laser  output  intensities  for  both  the  14A  and  21A  cases 
are  shown  in  Fig.  3:18.  The  output  shows  a gain  spike  at  the  onset  of 
the  laser  radiation  region  which  occurs  much  farther  downstream  in  the 
14A  case  due  to  the  lower  electrical  pumping  rate.  Of  course  the  gain 
spike  is  very  narrow,  thereby  containing  very  little  power,  and  will  only 
be  seen  in  an  oscillator  with  plane-parallel  mirrors.  The  radiation 
intensity  increases  beyond  this  gain  spike  until  the  electrical  excitation 
is  abruptly  ended  at  12.0  cm.  Combined  with  plane— parallel  sheet  mirrors 
which  are  as  high  and  as  long  as  the  electrical  excitation  region 
(8  x 11  x 140  cm),  these  results  predict  a power  output  of  45  kW  at  15.2% 
electrical  to  optical  output  efficiency  for  the  21A  case,  and  a power 
output  of  21  kW  at  11.4%  output  efficiency  for  the  14A  case. 

All  the  results  given  heretofore  are  for  a total  pressure  of 
280  Torr,  but  this  value  may  not  be  the  opLimum  pressure.  The  output 
efficiency  and  power  dependence  on  the  pressure  is  shown  in  Fig.  3:19 
for  the  14A  case.  The  efficiency  has  a peak  of  16.6%  at  90  Torr  and 
decreases  to  10.7%  at  300  Torr.  At  high  pressures  the  efficiency 
decreases  very  slowly  while  the  output  power  continues  to  increase  at  a 
large  rate.  Thus,  an  efficient  laser  using  the  1:7:20  gas  mixture 
should  operate  at  90  Torr  and  a high  power  laser  should  operate  at  high 
pressures . 

All  of  these  results  are  for  mirror  transmissivities  of  0%  and 
15%  with  absorption  and  scattering  losses  of  1.5%  at  each  mirror.  The 
output  power  and  efficiency  depend  on  the  mirror  transmissivity  and  this 
dependence  is  shown  in  Fig.  3:20  for  a current  of  21A.  These  results  are 
obtained  by  changing  only  the  transmissivity  of  the  15%  mirror.  A rather 
broad  peak  occurs  between  10%  and  15%  which  indicates  that  the  previous 
results  in  this  section  are  for  the  near  optimum  laser  output  coupling 
of  15%. 

A laser  should  operate  with  a radiaLion  intensity  within  the 
cavity  which  is  sufficiently  large  to  depopulate  the  upper  laser  level 
by  stimulated  emission  as  quickly  as  it  is  populated  (although  a steady- 
state  population  may  exist).  The  previous  results  for  vibrational 
temperatures  and  extractable  stored  energy  show  that  the  COFFEE  laser  is 
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Fig.  3:18  Calculated  laser  output  intensities  for  currents 
of  14  and  21  A. 
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indeed  operating  within  this  region.  Another  check  on  this  condition  is 
to  compare  the  intracavity  radiation  intensity  with  the  saturation 
radiation  intensity,  i . ihe  latter  intensity  is  that  intensity  which 
reduces  the  cavity  gain  to  one-hali  of  the  smali-sigial  gain.  By  assuming 
that  the  lower  laser  level  depopulation  is  independent  of  the  stimulated 
emission  rate,  1 is  that  intensity  at  which  the  deactivation  rate  of 

S cl  L 

the  upper  laser  level  by  stimulated  emission  is  exactly  equal  to  the 

98 

deactivation  rate  by  molecular  collisions.  The  stimulated  emission 

deactivation  rate  relative  to  the  collisional  deactivation  rate  is  given 

in  Fig.  3:21  as  a function  of  the  intracavity  intensity.  This  figure 

shows  that  for  the  given  conditions  the  saturation  intensity  for  the  COFFEE 

2 

laser  is  1.75  kW/cm  , which  is  indeed  much  smaller  than  the  calculated 

? 

operating  mtracavity  intensity  of  8.3  kW/cm  for  the  21A  case. 

The  examples  in  this  section  illustrate  the  utility  of  the 
CO 2 COFFEE  laser  kinetics  computer  code  and  demonstrate  the  excellent 
agreement  of  the  theoretical  small-signal  gains  with  the  experimental 
values.  This  code  can  be  used  either  to  explain  experimental  observations, 
or  to  develop  and  optimize  a laser  design  which  must  meet  specific 
application  requirements. 

3.1.5  Summary 

A harmonic  oscillator,  six-temperature  Boltzmann  equilibrium 
model  of  the  C0? : h'9 : He : H?0  COFFEE  laser  system  was  developed  and  tested. 

The  model  includes  the  gas  temperature  dependence  of  the  gas  flow  velocity, 
molecular  densities  (assuming  constant  pressure),  and  the  VV  and  VT 
collisional  relaxation  rates.  In  addition,  the  E/N  dependence  of  the 
electrical  excitation  rates  was  i.ncluded  as  well  as  the  effects  of  gas 
flow  and  external  mirrors.  All  important  sources  of  gas  heating  were 
incorporated,  including  direct  gas  discharge  heating,  VT  relaxation,  and 
intermode  W energy  defect  heating.  The  only  experimental  data  needed 
by  this  model  are  the  current  density  and  electric  field  distribution 
in  the  downstream  direction,  the  gas  mixture,  the  flow  velocity  and  gas 
temperature  at  the  entrance  of  the  laser  cavity,  and  the  geometrical 
design  of  the  laser  cavity. 
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The  model  is  quite  general  and  can  be  applied  to  most  C0?  gas 


flowing  laser  systems  of  present  interest.  However,  its  use  is  restricted 

to  electrical  excitation  and  stimulated  emission  rates  slower  than  the 

intraraode  relaxation  times  (^50  nsec),  to  systems  in  which  thermal 

diffusion  is  unimportant,  to  pressures  less  than  100  atm,  and  to  gas 

temperature  above  freezing  points  (^200°K)  but  below  dissociation  levels 

( u5000  K) . Furthermore,  the  basic  assumption  of  Boltzmann  equilibrium 

invalidates  the  model  at  high  radiation  intensities  or  very  low  CO 
29  2 

concentrations. 

Computer  predictions  were  compared  with  measurements  of  the 
P(22)  small-signal  gain  in  a 280  Torr,  1:7:20  flowing  gas  mixture  with 
2.0  Torr  H at  total  discharge  currents  of  1.0A  and  1.5A.  In  addition, 
these  same  conditions  were  modeled  with  external  resonators  included 
and  the  optical  cavity  theoretically  lengthened  to  140  cm,  although  no 
corresponding  experimental  energy  extraction  measurements  were  available 
for  comparison.  The  most  important  results  obtained  from  these  computer 
calculations  are: 

1.  The  model  gives  accurate  predictions.  The  six-temperature 
model  predicts  values  of  the  P(22)  small-signal  gain  which  agree  very  well 
with  the  values  measured  over  the  entire  cavity.  The  agreement  is 
particularly  good  near  the  peak  gain  where  the  experimental  measurements 
are  the  most  accurate. 

2.  Chemical  minority  species  are  important.  The  good 
agreement  between  theory  and  experiment  was  achieved  by  assuming  that 


the  laser  gas  contained  an  initial  background  concentration  of  0.15%  H20 
and  that  14%  of  the  was  chemically  converted  to  H 0. 


3.  The  gas  is  very  efficiently  pumped.  The  power  deposition 
into  the  gas  was  very  uniform  throughout  the  cavity,  and  the  operating 
value  of  E/N  was  always  close  to  the  optimum  value  for  efficient  pumping, 
vTith  75  to  83%  of  the  electrical  power  being  deposited  in  the 
vibrational  mode. 

4.  The  upper  laser  level  is  efficiently  depopulated.  The 
intracavity  radiation  intensity  was  about  five  times  larger  than  the 
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saturation  intensity  of  i.8  kW/cm  . This  assures  that  the  depopulation 
of  the  upper  laser  level  is  primarily  by  stimulated  emission  rather  than 
by  collisional  relaxation. 

5.  The  COFFEE  laser  is  efficient  and  powerful.  At  a pressure 
of  280  Torr  in  a 8 x 11  x 140  cm  cavity  an  ouLput  power  of  45  kW  is 
predicted  with  an  efficiency  of  15%.  Furthermore,  the  ouLput  efficiency 
peaks  at  a pressure  of  90  Torr,  whereas  as  the  output  power  peaks  at 

a pressure  above  300  Torr. 

6.  The  COFFEE  laser  can  be  easily  scaled.  The  predicted 

exit  temperature  for  the  8 x 11  x 140  cm  cavity  is  only  336°K  which  indicates 
that  the  laser  cavity  can  be  appreciably  extended  in  the  downstream 
direction  before  the  laser  performance  would  degrade  significantly  by 
collisional  deactivations  increased  by  temperatures  above  450°K. 

These  results  form  a firm  foundation  for  future  work.  The 
tests  of  the  computer  program  and  six-Lemperature  model  verify  their 
applicability  for  CO^  COFFEE  laser  kinetic  simulations  and  recommend 
their  use  for  system  modeling  and  prediction.  However,  detailed 
comparisons  with  experiments  under  lasing  conditions  are  required  for 
further  verification,  and  a better  understanding  of  the  hydrogen  to 
water  conversion  is  desirable.  With  these  reservations  the  CO^  COFFEE 
laser  kinetics  program  can  be  employed  for  systematic  variations  of 
parameters  to  establish  optimum  operating  conditions,  to  explain 
experimental  observations,  or  to  design  the  next  generation  COFFEE  laser, 

3.1.6  Conclusions 


Three  general  conclusions  can  be  drawn  from  the  laser  kinetics 
model  calculations: 

1.  The  six-temperature  model  can  accurately  predict  and  explain 
the  performance  of  CO^  COFFEE  lasers. 

2.  Efficient  (15%),  high-power  (45  kW)  COFFEE  lasers  can  be 
built  by  simply  increasing  the  present  optical  cavity  length  to  140  cm. 

3.  Appreciably  more  efficient  and  powerful  COFFEE  lasers  can 
be  built  by  extending  the  excitation  region  in  the  downstream  direction. 
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3.2  V-I  CHARACTERISTICS 


The  COFFEE  laser  excitaLion  scheme  uses  a continuous  self- 


sustained  glow  discharge  in  high-pressure  CO  laser  gas  mixtures.  The 


voltage-current  characteristics  of  the  glow  govern  the  power  into  the 
discharge  while  the  operating  E/N  controls  Lhe  laser  excitation  efficiency 
as  described  in  Section  3.1.  It  is  important,  therefore,  to  determine 
what  factors  govern  these  parameters. 

The  operating  E/N  of  a glow  discharge  is  Lhat  value  which 
leads  to  a steady  state  balance  between  the  rates  of  all  the  electron 


generation  and  loss  mechanisms.  The  analysis  of  Lowke,  Phelps,  and 
Irwin  shows  that  for  a volume  controlled  discharge 


j/N  = (eW  /y)(a-a)/N 


(3:9) 


where  j is  the  current  density,  N is  the  neutral  gas  density,  e is  the 
electron  charge,  W is  the  electron  drift  velocity,  y is  the  recombination 
loss  rate,  a is  the  ionization  coefficient,  and  a is  the  attachment 
coefficient.  Equation  (3:9)  expresses  j/N  as  a function  of  E/N  through 
the  dependencies  calculated  by  Lowke  et  al.  of  W,  a/N,  and  a/N  on  E/N. 

The  curves  in  Fig.  3:22  were  drawn  for  Eq.  (3:9)  using  three  representative 
gas  mixtures.  For  low  values  of  j/N  two  limits  are  shown.  The  broken 
curves  correspond  to  the  case  of  an  attachment  dominated  glow  while  the 
solid  curves  correspond  to  a recombination  controlled  glow.  Gases  such 
as  C02,  H20  and  02  are  known  to  produce  electron  loss  by  attachment.6 
For  example,  dissociative  attachment  by  CO  to  form  CO  and  0~  leads  to 


electron  loss.  Certain  gas  additives57-59  may  act  to  suppress  attachment 


and  lower  the  operating  E/N  at  a given  j/N.  This  can  result  in  higher 
laser  excitation  efficiencies.  The  operating  E/N  will  become  recombination- 
controlled  at  high  current  densities  since  recombination  losses  increase 
as  the  square  of  the  current. 
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The  effect  of  varying  the  helium  content  in  the  laser  mixture 
is  shown  in  Fig.  3:23.  If  the  proportion  of  helium  is  increased  for  a 
given  E/N,  the  average  electron  energy  is  also  increased  because  helium 
has  none  of  the  inelastic  vibrational  losses  that  are  present  in 
nitrogen  and  carbon  dioxide.  Thus,  electron  energies  sufficient  to 
excite  the  upper  laser  level  occur  at  lower  E/N  for  the  1:7:30  mixture 
than  for  the  1:7:0  mixture  where  helium  is  absent.  At  first  glance  one 
might  try  to  optimize  the  mixture  for  maximum  efficiency  by  varying  the 
helium  content.  However,  it  is  seen  from  Fig.  3:23  that  the  change  in 
E/N  for  maximum  efficiency  and  the  change  in  operating  E/N  are 
approximately  the  same.  It  appears,  therefore,  that  the  addition  of 
helium  has  no  significant  influence  on  the  effective  laser  excitation 
efficiency  for  a self-sustained  glow  discharge  of  interest. 

The  COFFEE  laser  operates  in  the  current  regime  (J /N  < 10~19 
A-cm)  where  attachment  appears  to  dominate  the  discharge.  Figure  3:24 
shows  calculated  values  of  a/N  and  a/N  as  a function  of  E/N  for  three  gas 
mixtures. ^ An  attachment  dominated  self-sustained  discharge  will 
operate  at  the  E/N  value  where  a = a.  Similar  calculations  for  a 1:7:20 
mixture  of  CO^N^He  predict  an  operating  E/N  of  3.0  x 10~16  V-cm2, 
which  is  independent  of  achievable  operating  current  densities.  This 
agrees  to  first  order  with  the  experimental  results  obtained. 

Ihere  are,  however,  a number  of  discrepancies  which  this 
simple  model  does  not  explain.  Figure  3:25  shows  the  predicted  and 
derived  E/N  distributions  as  a function  of  downstream  distance  in  the 
discharge  for  a typical  set  of  operating  conditions.  For  an  attachment 
dominated  discharge  in  a given  gas  composition  one  would  expect  the 
operating  E/N  to  be  a constant  independent  of  downstream  distance.  The 
derived  E/N  values  were  calculated  from  the  measured  E values  using  the 
laser  kinetics  code  to  determine  the  N distribution  as  described  in 
Section  3.1.  Near  the  upstream  portion  of  the  discharge  the  agreement 
was  within  10%.  However,  the  operating  E/N  dropped  almost  40%  as  the 
laser  gas  flowed  downstream  through  the  discharge,  which  is  presently 
unexplained  by  the  model.  If  the  discharge  were  recombination  controlled, 
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E/N  would  be  expected  to  increase  in  the  downstream  direction  because 
the  current  density  rose  (see  Fig.  3:3).  It  is  not  possible  to  explain 
a variation  of  this  magnitude  on  the  basis  that  the  discharge  may  be 
bending  more  at  an  upstream  location  than  at  a downstream  location. 

It  appears  that  something  is  created  in  the  upstream  portion  of  the 
discharge  which  is  carried  downstream  by  the  flow  and  effectively  lowers 
the  operating  E/N.  For  example,  various  transient  minor  molecular  species 
may  be  created  by  the  upstream  portion  of  the  discharge  which  subsequently 
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Such  processes 


alter  the  attachment-detachment  rates  downstream, 
have  not  been  modeled  yet. 

Two  other  parametric  dependencies  have  been  observed  for  E/N, 

“ Sh°™  ln  3:26  and  3:27>  «M<*  «•  not  fully  explained  at  present. 

Ihe  data  in  these  two  figures  are  typical  of  results  obtained  under  many 
operating  conditions  with  two  different  experimental  devices.  ^ The 
E/N  values  are  from  the  upstream  portion  of  the  discharge  where  gas 
heating  corrections  to  N did  not  have  to  be  taken  into  account. 

Figure  3:26  shows  that  Z/N  has  a slight  discharge  current  dependence. 
However,  the  E/N  dependence  on  gas  flow  velocity  shown  in  the  two  figures 
is  surprisingly  large.  For  example,  doubling  the  gas  velocity  increased 
E/N  by  v50,7.  Again,  this  cannot  be  explained  quantitatively  by  discharge 
del  lection  models.  Obviously  more  work  is  needed  in  this  area. 

In  summary,  it  should  be  noted  that  the  relatively  simple 
glow  discharge  model  does  give  a good  first  oruer  prediction  of  the 
operating  E/N.  In  addition,  the  operating  range  of  E/N  (Fig.  3:25) 

coincides  with  the  optimum  values  required  for  efficient  laser  excitation 
(Fig.  3:3). 
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4. 


H 1 GH  TOWEK  l.ASER  DEVELOPMENT 

One  oi  the  main  objectives  of  this  study  was  to  demonstrate 
the  feasibility  of  developing  multi-kilowatt  devices  using  the  COFFEE 
laser  excitation  technique.  The  discharge  scaling  experiments  described 
in  Section  2 and  the  modeling  studies  described  in  Section  3 clearly 
demonstrate  that  it  is  feasible  to  build  high-power  COFFEE  lasers.  In 
anticipation  of  this  result,  a parallel  effort  to  construct  a 25  kW 
proot  of  scalability  test  laser  was  started.  During  this  reporting 
period  a new  compact  closed-cycle  flow  facility  was  built  as  shown  in 
Fig.  4:1.  The  tank  is  6 feet  in  diameter  and  approximately  8 feet  long, 
i 1 w rates  up  to  <000  i/sec  are  attainable  with  operating  pressures  from 
50  to  760  Torr.  Heat  exchangers  and  an  external  DC  power  supply  (in 
back  room  of  fig.  4:1)  are  also  available  along  with  the  control  console 
shown. 

Preliminary  tests  using  an  existing  electrode  assembly  were 
performed.  The  maximum  usable  aperture  for  this  assembly  was  2 x 5 cm 
with  a 140  cm  active  path  length  along  the  optical  axis.  Up  to  6 kW 
*'f  output  power  was  obtained  from  this  discharge  volume  using  hole— 

' copied  optics.  The  conversion  efficiency  was  'vlO”'  with  ^4  W/cm  ^ being  extracted 
from  the  discharge  volume.  Calculations  using  the  kinetics  code  described 
in  Section  3 indicate  that  the  extraction  efficiency  of  the  simple 
optical  cavity  was  only  70’/. 

Future  plans  include  the  installation  of  a new  electrode 
assembly  with  a 10  x 10  cm  aperture  in  the  flow  system.  An  unstable 
resonator  optical  cavity  will  also  be  installed.  Originally  this  work 
was  to  be  performed  as  a Phase  III  continuation  of  Phase  II  work  described 
in  this  report.  However,  the  lack  of  funding  has  for  the  time  being, 
delaved  this  continuation. 
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Fig.  4:1  Compact  COFFEE  laser  facility  for  25  kW  scaling  tests 


5 . SUGGESTED  FUTURE  WORK 


r. 

n 
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The  scaling  criteria  for  the  COFFEE  laser  excitation  scheme 
are  developed  to  the  point  that  simple,  rugged,  and  compact  high-power 
laser  systems  can  easily  be  designed.  There  are  a few  areas,  however, 
where  additional  information  would  be  most  useful.  First,  certain  portions 
of  the  discharge  physics  need  to  be  developed  more  fully.  For  exampie, 
the  effects  of  transient  minor  molecular  species  and  gas  flow,  as  discussed 
in  Section  3.2,  on  the  operating  E/N  are  poorly  understood.  Additional 
knowledge  in  this  area  would  better  .v  ire  controlled  operation  at  the 
optimum  E/N  for  efficient  excitation.  . urthermore,  the  glow-to-arc 
transition  problems  appear  to  be  affected  by  the  same  processes.  Any 
information  which  helps  raise  the  arcing  limit  will  result  in  increased 
specific  input  energy,  input  power  density,  efficiency,  and  power  output. 

Second,  a large  scale  device  should  be  built  and  tested,  not 
only  to  assess  the  reliability  of  the  scaling  criteria,  but  also  to 
determine  if  there  are  any  major  engineering  problems  which  have  to  be 
overcome.  Although  no  significant  problems  are  anticipated,  proof  of 
scalability  should  be  demonstrated. 

Finally,  the  cavity  extraction  and  beam  quality  problems, 
which  tend  to  be  inherent  in  all  high-power  cw  devices,  should  be  addressed 
more  fully.  The  COFFEE  laser  has  the  unique  feature  that  the  input  power 
is  easily  tailored  in  the  transverse  direction  in  order  to  match  certain 
cavity  extraction  criteria.  All  this  additional  information  will  provide 
a better  technology  base  for  designing  the  next  generation  devices. 
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7.  APPENDICES 


APPENDIX  A 

NOMENCLATURE  FOR  SECTION  3.1 


Boltzmann  Factors 


T 

= 

gas  kinetic 

T 

S 

vibration.il 
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e-VTb 
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= 

u-0n/T„ 
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-0- /Ta 
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u---u  lw 
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c-->b/T 

n 

sc 
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= 

e--w/T 

temperature 

temperature  of  u vibrational  mode 


Characteristic  Temperatures 

= characteristic  temperature  of  a vibrational  mode 
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°a 
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3380°K 

*n 
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3354° 

°w 
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2295° 
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= 

1997° 

°b 

= 

960° 

wb 
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6w 
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0b  = 1335 

aw 
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0 = 1085 
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nw 

= 

9 
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- 

<j  = 1059 
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ab 
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0 
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36L  = 599 
b 

sb 

= 

Gs 

- 

2>j  = 147 

b 

an 

= 

Ja 

- 

•)  = 26° 
n 

= characteristic  temperature  of  rotational  mode 


0. 556°K. 
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Energy  Defect  Factors 
A = e"°an/1 
B - e-^b/T 
C = e“Jab/T 
n . „-0nw/T 


D = e nw/ 1 
F = e“'Jaw/1 


G = e 


-,Jwo/T 


Gas  Heating  Factors 


Cy  = specific  heat  at  constant  pressure  in  joule  deg  * mole~\ 

Cu  = 20.79 
h 

C = 29.11 
n 

C = 29.11 
c 

C = 33.26 

w 

E^  = characteristic  energy  density  in  joule  cm"3. 

Fe  “ fraction  of  the  electrical  energy  which  is  not  lost  by 
photon  emission  from  electronically  excited  levels, 
fa  = molar  density  in  no.  moles  cm  3 of  a species. 

C = electronic  heating  factor— the  fraction  of  the  electronic 
excitation  energy  which  heats  the  gas. 


Radiation  Intensities 


aiYj  = loss  Per  bounce  at  the  ith  mirror  (including  transmis- 
sivities) . 

byj  = scattering  loss  per  cm  inside  the  gas. 

= total  light  intensity  in  warts/cm  . 

= collimated  light  intensity. 

1^1  = uncollimated  light  intensity. 

IQ  = incident  light  intensity  lo  be  amplified  by  laser. 

Iout  = total  light  intensity  out  of  the  laser. 

K = V^/A^  = uncollimated  light  conversion  factor  from  power 
density  to  light  intensity  in  cm. 

K"  = Ao/Vc  = input  1:Lght  conversion  factor  from  light  intensity 
to  power  density  in  cm-3. 
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PyJ  = sponCaneous  emission  power  density  in  watts/cmJ. 
ayj  total  gain  of  gas  medium  in  cm  \ 

l\j  = ayerage  cavity  loss  per  cm  (also  called  the  cavity 
gain  or  threshold  gain). 
oI\j  = P°Pu-lation  inversion  density  in  cm~3. 

£yj  ” i ratal i ] i ty  that  a spontaneously  generated  photon  is 
led  into  the  Jasing  mode. 

‘iyj  “ transmissivity  of  the  ith  mirror. 

Subscripts 

a = asymmetric  stretch  mode  (ASM)  of  CO^. 
b = bending  mode  t,BM)  of  CO^. 

c * CO,,, 

h = He. 

n = N . 

o = denotes  value  for  the  gas  flowing  into  the  discharge 

volume  before  that  value  is  modified  by  the  discharge 
or  lasing  action,  i.e.,  the  value  at  x = 0. 
s = symmetric  stretch  mode  (SSM)  of  CO  . 
w = H^O. 

Transition  Rates 
Collisional 

A\B\C'  = parameterization  constants  for  reaction  rate  constants 

^ = = constant  used  in  reaction  rate  constant. 

kaB  " reaction  rate  constant  in  cm^  per  3 particle  per  second. 

■ no.  of  a particles  cm  ^ (Note  that  N =N  =N  =N  ) 

a b s c 

N = total  number  of  density  of  gas  in  particles  cm-3. 

Ra  = VT  transition  rate  in  no.  of  transitions  per  second. 

Ra0  ~ W transition  rate  due  to  collisions  between  a and  $ 
particles. 
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Electrical 

E ” electric  field  in  V'  in. 

I = discharge  current  in  A. 
j = discharge  current  density  in  A cm-". 

= electrical  pumping  term  ir.  sec 
V = voltage  drop  across  the  discharge  in  volts. 

W = power  density  entering  the  a mode  in  W cm-3. 

Y0  ~ fraction  of  electrical  power  density  which  excites  the 
i vibrational  mode,  electronic  levels  (a=e) , or  ionizes 
the  gis  (a=i). 


Stimulated  Emission 
A 


21 


YJ 

fYJ 


= spontaneous  transition  probability  rate  between  upper 
and  lower  laser  levels  in  sec 

= Gaussian  bandwidth  in  cm  1 (full  linewidth  at  half  maximum). 

= Lorentzian  bandwidth  in  cm  1 (full  linewidth  at  half  maximum) 

— 1 —2 

= photon  flux  in  photons  sec  cm 

= rotational  quantum  number  of  upper  laser  level  (must  be 
odd) . 

= rotational  quantum  number  of  lower  laser  level  (must  be 
even)  . 


CO^  mass  in  amu. 


Yd 

"Yd 


, - population  density  of  upper  laser  level  in  particles  cm 
population  density  of  lower  laser  level  in  particles  cm 


Yd 

"ci 


-3 

i 

-3 


Yd 


syJ 

21 


= denotes  P branch  transition  (J  = J-l) . 

upper 

= common  rotational  factor  of  the  stimulated  emission  factor. 

“ denotes  R branch  transition  (J  = j+1). 

upper 

- stimulated  transition  probability  rate  in  sec-1. 

= fraction  of  molecules  in  the  Jth  rotational  level. 

= P or  R 

= stimulated  emission  term  in  sec 
= stimulated  emission  factor 


specifies  branch  for  rotational  transition. 

-1 


reduced  mass  of  species  C02  and  i in  amu. 
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v r-  laser  radiation  wave  number  in  on  1 . 
o 

,(0)  = value  of  line  shape  factor  't  line  center  in  cm. 


yJ 


Cu 


yJ 


= stimulated  emission  cross-section  in  cm^ 
= optical  broadening  cross-section. 

o 

o = 91 
cn  o 

o = 25  A 
ch  o 

a = 130  A 

CC  o 

o = 73  A 

ew 

= Lorentzian  lifetime  in  seconds. 


Laser  Cavity 

The  following  relations  (in  cgs  units)  pertain  to  the 
geometrical  design  of  the  laser  cavity  and  are  easily  derivable.  A 
plane-parallel  laser  cavity  with  a rectangular  discharge  is  assumed. 

The  mirrors  are  assumed  to  be  rectangular  and  plane-parallel. 

x = downstream  distance  from  the  entrance  to  the  laser  cavity, 
v - gas  flow  velocity  in  cm  sec*1. 

L • distance  between  the  mirrors. 


dis 


l 

W 


* WH1  = volume  of  the  gas  discharge. 

= length  of  discharge  along  a line  joining  the  mirror  centers. 

= length  of  discharge  perpendicular  to  the  line  joining 
the  mirror  centers  and  along  the  direction  of  the  gas 
flow. 

= length  of  discharge  between  the  electrodes,  and  the 
height  of  the  mirrors. 

cross-sectional  area  of  input  radiation  intensity. 

■“  = total  volume  from  which  the  uncollimated 

radiation  originates. 

= 2£(W+H)  = enclosing  area  of  discharge  volume  through 
which  the  uncollimated  radiation  escapes. 

= Vdis  “ t0  tal  volume  of  gas  which  contains  all  the 
lasing  radiation. 
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Cons  tants 


k = 1.38054 

k'  = 1.38054 

h'  = 6.6256  x 

c - 2.997925 

e = 1.60210  : 

tt  = 3.141593 

N = 6.02252  ; 
o 

W - 1.65979  ) 


. ln-16  , -l 

• 10  erg  deg 

: 10  ^ joule  deg 
-34 

10  joule  sec 

m10  “I 

x 10  era  sec 

-19 

10  coulombs 


10 

10 


23 

-24 


particles  mole 
-1 

gm  amu 
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APPENDIX  B 


COLLISIONAL  RELAXATION 


The  collisional  exchange  of  energy  between  molecules  in 
different  vibrational  states  has  been  treated  by  the  SSH  theory, 20 
which  employs  the  Landau-Teller  selection  rules^^  for  quantum  transitions 
of  simple  harmonic  oscillators.  The  dominant  collision  si  interactions 
are  depicted  in  Fig.  3:2.  The  A,  B,  C,  D,  F and  G factors  in  the  VV 
relaxation  terms  of  Eqs.  (3:1-6)  account  for  the  energy  defects 
involved  in  the  collisional  interactions. 

The  double  degeneracy  of  the  CO^  bending  mode  gives  rise  to 
the  first  and  last  two  factors  o^  2 in  Eq.  (3:3).  In  keeping  with  the 
SSH  theory  approximations,  only  those  collisional  transitions  are  assumed 
to  occur  which  involve  one  of  the  bending  modes,  independent  of  the  other 
bending  mode.  The  factor  of  3 arises  from  the  exchange  of  one  ASM  quantun 
for  three  BM  quanta,  while  the  second  factor  of  2 arises  from  the  exchange 
of  one  SSM  quantum  for  two  BM  quanta.  All  other  VV  energy  exchanges 
involve  one  quantum  in  each  mode. 


The  total  VT 

transition 

*b 

= E 

K 

k6 * (8) 

kb 

R 

P 

= E 

(3) 

k 

n 

R 

3 

= E 

n 

(3) 

k 

w 

3 

w 

3 = c,h,n,w 


(B : 9) 


(6)  3 

where  is  the  reaction  rate  constant  in  cm"5  per  particle  3 per 

second,  and  N^_  is  the  number  density  of  the  particles.  The  sum  is 

over  all  the  gas  components  S(c=CO,,,  h=He,  n=N9,  w=H90) . The  W 

transition  rates  for  intramolecular  energy  exchanges  through  collisions 

with  other  molecules  are  given  by 
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(B: 10) 


R , = E 
ab  „ ab 

3 

R = 1 k^  N 
sb  sb 


whereas  the  intermolecular  VV  transition  rates  are  given  by 


an 


k N 
an  n 


na 


k N 
na  a 


aw 

^w 

R 

nw 


k N 
aw  w 


1c  N 

DW  W 


k N 
nw  w 
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k N 
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wb  b 


k N 
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(B : 11) 


In  these  expressions 


k = k 

ub  3a 


(B:12) 


and  the  a,  b,  and  s subscripts  denote  the  ASM,  BM,  and  SSM  of  the  CO, 

l- 

molecule.  The  particle  densities  are  speci  ied  by  the  laser  operating 
conditions,  and  the  reaction  rate  constants  are  derived  from  experimental 
measurements  or  estimated. 

Qualitatively,  collisional  relaxation  rate  constants  are 
fairly  well  understood  both  theoretically  and  experimentally.  However, 
quantitatively  the  theory  is  inadequate  and  the  experiments  aic  frequently 
very  difficult  to  perform.  .Nonetheless,  many  relevant  rate  constants 
have  been  measured.  For  these  reasons  the  theory  is  used  to  give  a 
parameterized  equation  for  the  temperature  dependence  which  is  then  fitted 
to  available  measurements. 

30-32 

Collisional  theory  states  that  the  parameterized 

temperature  behavior  of  W and  VT  collisional  rate  constants  is  given  by 


A'x2  -B'x+C'x^D'x2 
3 e 

1 + F'xJ 


(B:  13) 
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where 


x = 


t-1/3 


and  the  A',  B',  C',  D',  r- ' F*  are  parameterization  constants.  Actually, 
D'  is  fairly  well  determined  to  be 


where  0 is  the  energy  defect  of  the  exothermic  collisional  transition. 
Since  the  temperature  dependence  of  Eq.  ( B:  13)  is  dominated  by  the  ter^Q 
in  x in  the  exponent,  a reasonable  simplification  is  to  take 


F'  = 0 


(B:  15) 


Equation  (B:13)  was  fitted  to  experimentally  measured 
values  at  three  different  temperatures  by  using  Eqs.  ( B: 14)  and  (B:15), 
This  fit  then  determined  the  values  of  A',  B'  and  C'  and  the  collisional 
rate  constants  used  in  the  laser  kinetics  code. 

This  fitting  procedure  can  be  used  in  two  different  ways.  If 
B'  and  C’  are  positive  and  have  reasonable  values,  then  the  theory  uoon 
which  Eq.  (B:13)  is  based  is  probably  valid.  In  this  case 


A' 

B' 

C' 


7/6 


1/3 


1/2  1/6 

e U 


(B : 16) 
(B: 17) 
(B: 18) 


where  p is  the  reduced  mass  of  the  colliding  particles  and  e is  the 
Lennard-Jones  intermolecular  potential  energy  parameter.  ^ These 
equations  can  then  be  used  to  scale  known  rate  constants  for  an  estimate 
of  unknown  rate  constants.  On  the  other  hand,  if  B'  or  C'  is  negative 
or  unreasonably  valued,  Eq.  (B:13)  can  be  used  simply  as  a mathematical 
fitting  expression,  and  scaling  is  invalid. 
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A criticaJ  ^evaluation  of  available  rate  constants  shows  that 
Iaylor  and  Bitterman  have  <i.,ojmbled  the  most  reliable  set  of  experimental 
measurements.  All  relevant  VT  rate  constant  parameters  can  be  derived 
from  this  reference  and  are  given  in  Table  B:l. 

The  VT  rate  constants  for  the  relaxation  of  the  BM  are  shown 
in  Fig.  B: 2 2.  Water  clearly  exhibits  anomalous  behavior  in  that  its  rate 
constant  for  emptying  the  lower  laser  level  (since  the  BM  and  SSM  are 
tightly  coupled)  is  highest  at  low  temperatures. 

The  VT  rate  constants  for  the  relaxation  of  the  NVM  are  shown 
in  Fig.  B: 23.  The  relaxation  by  C02  was  scaled  (see  Eqs.  (B:16-18))  from 
the  experimental  results  for  relaxation  by  N . Water  is  clearly  the 
most  effective  relaxant  of  Since  this  relaxation  competes  with  the 

W pumping  of  CO 2 by  excited  N^,  high  concentrations  of  H90  should  be 
avoided  when  possible.  On  the  other  hand,  1^0  itself  can  pump  C02  by 
VV  collisional  interactions,  although  the  pumping  is  less  efficient. 

Clearly  the  t^O  concentration  plays  an  important  role  in  C02  laser  kinetics, 
particularly  in  the  lower  temperature  range. 

The  only  measured  VT  rate  constant  for  the  relaxation  of  the 
WBM  is  that  due  to  other  H20  molecules.  The  result  is  plotted  in  Fig. 

B:  2^*  However,  the  data  points  (Fig.  10  in  Ref.  25  ) separate  into  two 
groups  which  have  such  a large  spread  that  even  the  sign  of  the  slope 
at  500  K is  uncertain.  The  values  themselves  should  not  be  trusted 
outside  the  temperature  range  of  300  to  600°K.  The  other  rate  constants 
have  been  estimated  by  scaling,  and  by  arbitrarily  reducing  the  peak  values 
to  provide  consistency  with  the  abnormally  rapid  deactivation  by  II20 
itself.  This  procedure  can  yield  only  order-of-magnitude  estimates  at 
best,  and  is  employed  as  a matter  of  necessity  until  better  experimental 
values  are  available.  In  any  event,  the  low  concentration  of  H20  in  laser 
mixtures  reduces  the  importance  of  accurate  rate  constants  for  H20 
deactivation  processes. 

1 
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TABLE  B.  1 


VT  Rate  Constant  Parameters  for  k 

i 


i 

a 

A' 

B’ 

C' 

D' 

Source* 

b 

C°2 

4.99 

(-05) 

191 

348 

480 

TB 

N2 

9.68 

(-06) 

192 

358 

480 

TB 

He 

2.67 

(-06) 

121 

149 

480 

TB 

h2o 

7.93 

(-12) 

- 58 

-267 

.80 

TB 

n 

C°2 

1.09 

(-03) 

276 

- 42 

1677 

ES 

N2 

8.03 

(-04) 

254 

-646 

1677 

TB 

He 

4.44 

(-07) 

116 

-350 

1677 

TB 

H2° 

6.56 

(-05) 

123 

-386 

1677 

TB 

w 

C°2 

1.24 

(-19) 

-308 

-1446 

1147 

ES 

N2 

1.19 

(-17) 

-315 

-1856 

1147 

ES 

He 

5.18 

(-18) 

-240 

-1111 

1147 

ES 

h2o 

7.59 

(-15) 

-268 

-1468 

1147 

TB 

TB  - Taylor  and  Bitterman  ^ 
ES  = Estimated 
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25 

Taylor  and  Bitterman  have  also  compiled  the  best  set  of 
experimental  measurements  for  W rate  constants.  The  relevant  W 
rate  constant  parameters  are  given  in  Table  B.2.  Unfortunately,  the 
W rate  constant  in  which  helium  is  the  collision  partner  lias  not  been 

•J  O 

measured,  but  it  can  be  estimated  by  scaling  the  results  of  MacDonald 
for  argon.  These  rate  constants  are  plotted  in  Fig.  B:25  where  the 
anomalous  behavior  of  H90  is  evident  once  again.  The  high  values  of  H20 
indicate  that  the  H^O  concentration  should  be  kept  low  since  it  efficiently 


depopulates  the  upper  laser  level. 

Measurements  of  the  VV  relaxation  rate  constants  between  the 

S'JM  and  BM  have  been  the  source  of  considerable  controversy  since  implied 

rates  of  emptying  the  lower  laser  level  differ  by  more  than  an  order  of 

magnitude.  ‘ The  values  implied  by  Bulthuis  and  Ponsen^  were  at  500°K 

which  precludes  the  deduction  of  any  temperature  dependence.  However, 

the  main  problem  is  that  the  values  are  thirty  times  slower  than  tnose 

deduced  by  Sharma  from  the  measurements  of  Rhodes,  Kelly,  and  Javan. 

More  recent  experiments  lave  been  easier  to  interpret  and  have  veriiied 

37 

the  more  rapid  rate  consent.  In  particular,  Stark  has  measured  the 
rapid  relaxation  by  C0,(,  and  He  at  400°K.  Although  Stark  gives  no 


temperature  dependence,  these  values  can  be  combined  with  the  temperature 

38  r 

dependence  measured  by  Bulthuis  for  the  relaxation  by  C07  (the  Bulthuis 

measurements  give  values  which  are  too  low  to  use  directly).  The  relaxation 

by  HO  is  estimated  by  again  employing  Bulthius’s  temperature  dependence 

39 

with  the  measurement  of  Rosser  et  al.  The  results  of  these  estimates 
are  shown  in  Fig.  B:26.  The  collisional  relaxation  between  the  SSM  and 
BM  is  obviously  much  faster  than  that  between  other  modes,  as  expected 
from  the  near  resonance  of  the  two  modes. 


The  intermolecular  W relaxation  rate  constants  are  shown  in 

25 

Fig.  B:27.  The  k values  are  taken  from  Taylor  and  Bitterman,  whereas 
an  25 

the  values  are  taken  from  Taylor  and  Bitterman  at  high  temperatures 
and  from  Evans  at  low  temperatures.  The  values  for  k are  estimated  bv 
scaling  k^  by  the  reduced  masses  only.  Attempts  to  measure  have  been 
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TABLE  B.2 


W Rate 

Constant  Parameters 

tor  kjj 

(a)  and 

kij 

i 

1 

a 

A’ 

B* 

C' 

D' 

Source* 

ASM 

BM 

co2 

1.45  (-02) 

298 

830 

300 

TB 

N2 

5.12  (-02) 

330 

952 

300 

TB 

He 

3.24  (-02) 

330 

952 

300 

ES 

h20 

2.50  (-12) 

- 48 

-270 

300 

TB 

SSM 

BM 

C°2 

3.26  (00) 

304 

974 

74 

ES 

N2 

1.50  (00) 

304 

974 

74 

ES 

He 

1.76  (00) 

304 

974 

74 

ES 

HO 

2.30  (00) 

304 

974 

74 

ES 

*4* 

ASM 

NVM 

— 

4.07  (-10) 

368 

115 

13 

TB 

ASM 

NVM** 

— 

1.19  (03) 

602 

2895 

13 

TB 

ASM 

WBM 

— 

8.00  (-04) 

215 

430 

542 

ES 

BM 

WBM 

— 

7.93  (-14) 

- 58 

-267 

668 

ES 

NVM 

WMB 

— 

9.57  (-04) 

226 

442 

530 

TB 

* 25 

TB  * Taylor  ar.d  BitLennan 

ES  » Estimated 

+Valid  below  1000°K 

**Valid  above  1000°K 
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unsuccessful,  probably  because  this  process  Is  masked  by  the  k^W^  VI 
process;  was  arbitrarily  assured  to  be  0.G1  k^w) . The  results  In 
Fig.  B:27  Indicate  that  water  is  desirable  at  low  temperatures  but 
undesirable  at  high  temperatures,  and  that  nitrogen  pumps  the  carbon 
dioxide  more  efficiently  et  low  and  high  temperatures  than  it  does  at 
intermediate  temperatures. 

The  results  of  this  section  clearly  show  that  much  experimental 
work  remains  to  be  done  in  the  determination  of  VT  and  W rate  constants. 
Data  is  especially  lacking  on  all  interactions  which  involve  water. 

This  deficiency  is  unfortunate  because  the  temperature  dependence  of 
water  rates  is  known  to  be  anomalous  and  the  rates  abnormally  high,  and 
hence  unscalable  with  present  theoretical  models.  Moreover,  small 
concentrations  of  water  are  known  to  have  significant  effects  upon  the 
behavior  of  CC>2  laser  devices. 
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APPENDIX  C 


\ I 


. 3 


ELECTRICAL  EXCITATION 

The  electrical  pumping  rates  are  easily  determined  from 
Fig.  3:3.  1 he  total  electrical  power  density  into  the  discharge  is 

given  by  the  product  of  the  discharge  current  density  j and  the 
electric  field  E across  the  discharge  as  measured  experimentally.^ 

The  electrical  power  density  entering  the  a vibrational  mode  is  given  by 


W = y j£ 

u or 


(C:19) 


where  is  the  fraction  ol  electrical  power  density  which  enters  the 
a vibrational  mode  (plotted  as  a percentage  in  Fig.  3:3).  The  transition 
rate  into  a particular  mode  is  given  by  the  power  density  input  W , 
divided  by  the  energy  level  spacing  k'Q  and  the  number  of  a molecules 
per  unit  volume  traveling  in  the  flow  direction: 


u = w /k'e  n 

a a a a 


( C:  20) 


evaluated  at  the  location  x.  This  expression  is  the  electrical  pumping 
term  contained  in  the  basic  rate  Eqs.  (3:1-6),  and  is  valid  provided  that 
(1)  the  change  in  the  discharge  parameters  as  the  gas  flows  downstream 
occurs  at  a slower  rate  than  the  electron-neutral  inelastic  collision 
rates,  (2)  at  any  given  downstream  location  the  discharge  is  spatially 
uniform  in  the  plane  defined  by  line  of  sights  between  the  mirrors  and 
between  the  electrodes,  i.e.,  the  plane  perpendicular  to  the  direction 
oi  rlow,  and  (3)  all  molecules  are  in  the  ground  state  of  internal  energy.  ^ 
These  conditions  are  satisfied  in  most  CO,,  lasers. 


114 


A quantity  of  interest  related  to  the  electrical  excitation  is 
the  maximum  extractable  stored  energy  density  in  the  laser,  which  is 

_3 

defined  in  joule  cm  as 


E = (~ ) E + (— ) E 
s v ' a v n 
a n 


(C:21) 


E = N k ' 0 
a a a 


( C: 22) 


where  6 is  the  characteristic  temperature  of  the  a mode,  E and  E ar° 
u a n 

the  characteristic  energy  densities  for  the  CO^  asymmetric  mode  and  the 

n2  vibrational  mode  respectively,  and  the  v's  are  the  appropriate 

frequencies.  This  definition  is  realistic  if  the  equilibration  of  the 

set  of  rotational  levels  belonging  to  a gi/Ln  vibrational  level  is 

instantaneous  usually  a good  approximation.  Within  this  approximation, 

all  the  rotational  levels  associated  with  a vibrational  level  participate 

9 / 

in  a given  vibrational- rotational  transition.  Equation  (3:21)  can  be 
expressed  in  terms  of  the  Boltzmann  factors  as 


E = h'c  v — N + N ] 

o 1-a  c 1-n  n 


( C: 23) 


where  is  the  wave  number  of  the  laser  radiation.  This  expression  is 
the  total  energy  density  stored  in  C0?  asymmetric  stretch  mode  and  in  the 
vibrational  mode  which  can  be  extracted  at  10.6  _,m,  assuming  one  photon 
for  each  vibrational  quantum. 

Stored  energy  densities  in  C02  lasers  are  very  high  and  explain 
why  they  are  so  efficient  and  have  such  high  output  powers.  However,  stored 
energy  is  only  half  the  story— it  is  worthless  if  it  cannot  be  depleted  by 
stimulated  emission.  This  fact  explains  why  it  is  difficult  to  optimize 
laser  operation— a change  in  operating  parameters  may  increase  the  stored 
energy  density  but  decrease  the  stimulated  emission  depletion  of  that  energy, 
or  vice  versa. 
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APPENDIX  D 


GAS  HEATING 
42 

It  is  well  known  that  the  kinetic  temperature  of  the  gas  is 
important  in  molecular  lasers.  The  collisional  rates  change  as  the  gas 
temperature  varies  and  as  shown  in  Appendix  B these  rate  changes  can  be 
very  large.  Furthermore,  the  flow  velocity  and  gas  density  depend 
directly  on  the  gas  temperature.  It  is  therefore  imperative  to  determine 
the  kinetic  gas  temperature. 

Equation  (3:6)  is  derived  by  applying  the  ideal  gas  low  to  the 
specific  heat  equation  at  constant  pressure.  43  The  left-hand  side  gives 
the  heat  energy  per  unit  volume  absorbed  by  the  gas  and  is  determined  by 

the  specific  heats  of  the  gas  components,  each  component  a weighted  by  its 
molar  density: 


The^  specific  heats  of  the  gas  components  at  constant  pressure  are  given 
by 


Ch  - 20.79 

C = 29.11 
c 


C = 33,26 
w 

in  joules  deg  1 mole  1.  These  values  are  determined  from  the  full 
excitation  of  the  translational  and  rotational  degrees  of  freedom  of  the 
molecule  with  no  excitation  of  the  vibrational  degrees  of  freedom.  45 

This  is  an  excellent  approximation  in  the  gas  temperature  range  encountered 
in  CO^  gas  lasers. 
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The  specific  heats  at  constant  pressure  are  certainly  the  values 

to  use  for  the  COFFEE  laser,  and  are  probably  the  best  values  for  long 

pulse  lasers  where  the  gas  has  tine  to  expand,  and  for  high  energy  lasers 

where  the  rapid  electrical  heating  of  the  gas  causes  a rapid  expansion.  ^ 

However,  for  very  short  pulses  the  gas  does  not  have  time  to  expand  and 

23 

the  specific  heats  at  constant  volume  should  be  used. 

As  discussed  above,  the  left-hand  side  of  Eq . (3:6)  describes 
the  distribution  of  heat  energy  into  the  heat  content  of  the  molecules. 

The  molecular  kinetic  energy  is  described  by  the  gas  temperature.  The 
right-hand  side  describes  the  sources  and  sinks  of  the  heating.  There 
are  three  major  sources  of  heating:  energy  defect  heating  from  W 

collisional  transitions,  gas  heating  from  VT  collisional  transitions, 
and  direct  heating  from  the  gas  discharge.  All  three  sources  are  included 
in  the  present  model.  However,  heat  sinks  such  as  cold  walls  are  not 
included.  The  latter  approximation  is  valid  as  long  as  the  flow  velocity 
is  not  so  low  that  thermal  diffusion  to  the  enclosure  walls  becomes 
important. 

The  energy  defect  heating  terms  from  W collisional  transitions 
are  given  on  the  right-hand  side  of  the  first  three  lines  of  Eq.  (3:6). 

The  gas  heating  terms  from  VT  collisional  transitions  are  given  on  the 
fourth  line  of  Eq.  (3.6).  In  a typical  COFFEE  laser  discharge  W heating 
is  less  than  10%  of  the  total  heating  whereas  VT  heating  is  over  80%  of 
the  total  heating  (see  Fig.  3:12  and  3:13). 

The  electrical  energy  density  available  to  the  gas  is  jE  where 
the  electric  field  E and  current  density  j are  measured  quantities. 

However,  not  all  of  the  electrical  energy  is  available  to  heat  the  gas 
directly  a good  portion  of  the  electrical  energy  excites  molecular 
vibrations  rather  than  increasing  the  translational  energy  of  the 
molecules.  The  vibrational  excitation  energy  must  therefore  be  deducted 
from  the  total  electrical  energy.  Some  of  this  vibrational  excitation 
energy  will  eventually  heat  the  gas  through  the  VV  and  VT  relaxation 
processes  which  have  already  been  included.  In  the  last  line  of  Eq.  (3:6), 


the  first  term  (.without  the  factor  F^)  is  the  electrical  energy  density 
available  to  the  gas  and  the  terms  enclosed  in  parentheses  deduct 

the  electrical  energy  density  expended  in  vibrational  excitation  rather 
than  in  gas  heating. 

ihe  factor  F£  is  the  fraction  of  the  electrical  energy  density 
which  is  not  lost  to  the  system  through  the  electronic  excitation  of  the 
C°2  and  N2  molecules  which  then  relax  by  photon  emission,  rather  than  by 
relaxing  through  collisions  which  heat  the  gas.  Since  not  all  electronically 
excited  C02  and  N'2  molecules  relax  by  photon  emission,  the  factor  F can 
be  written  as 


Fe  " 1 “ (1~Oye  ( D:  26) 

where  is  the  fraction  of  the  electrical  energy  which  excites  the  gas 
into  electronic  levels  (see  Fig.  3:3)  and  t,  is  the  fraction  of  this 
excitation  energy  which  heats  the  gas  directly.  The  factor  ; is  called 
the  electronic  heating  factor,  and  1-;  gives  the  fraction  of  the  electronic 
excitation  energy  which  escapes  the  laser  system  by  high  energy  photon 
emission.  Since  the  electronic  heating  factor  is  not  included  in  other 
laser  kinetic  models,  it  will  be  examined  in  detail. 

The  electrical  energv  supplied  to  a CO^N^He  laser  mixture  is 
divided  among  elastic  and  rotational  processes,  vibrational  and  electronic 
level  excitations,  and  ionization.  The  exact  division  depends  on  the 
relative  concentrations  of  the  gases  and  the  operating  values  of  E/N.  For 
example,  in  a 1:7:20  laser  mixture  at  an  E/N  value  of  2 x 10~16  V-cm2 
(see  Fig.  3:3),  83%  of  the  electrical  energy  is  deposited  into  the  C02 
00  1 level  and  the  lowest  eight  vibrational  levels  of  N2>  and  the  remaining 
17  goes  into  elastic  losses,  rotational  excitation,  and  the  excitation 
of  the  C02  bending  and  symmetric  stretch  vibrations.^7  There  is  only  5%  of  the 
energy  in  electronic  excitation  or  ionization. 

At  higher  values  of  E/N  significant  percentages  of  the  electrical 

energy  will  excite  the  C02  and  N2  electronically.  As  shown  in  Fig.  3:3, 

this  percentage  is  near  zero  for  E/N  - 1.0  x 10“16  V-cm2  but  rises  rapidly 

to  82%  near  10  V-cm2.  For  CO  lasers  which  operate  at  E/N  values  below 
^ ^ ^ ^ ^ 

2.0  >.  10  V-cm  , it  is  clear  from  Fig.  3:3  that  electronic  excitation 
can  be  neglected  and  previous  laser  kinetics  models  are  quite  applicable 
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in  this  range,  which  includes  most  CO  lasers.  Non-flowing,  self-sustained 

^ — 16  2 

1:7:20  pulsed  discharges  normally  operate  at  3.0  x 10  V-cin  where 

electronic  excitation  cannot  be  neglecLed,  whereas  the  COFFEE  laser  operates 
“16  2 

around  2.0  x 10  V-cm  where  the  electronic  excitation  is  quite  small. 

Thus  the  electronic  heating  factor  cannot  be  neglected  in  pulsed  self- 
ustained  discharges  but  is  relatively  unimportant  in  the  COFFEE  laser. 

For  this  reason  an  accurate  value  of  the  electronic  heating  factor  is 
not  necessary  and  the  empirical  values  of  0.40  will  be  taken  since  this 
value  works  so  well  for  self-sustained  discharge  pulsed  CO^  lasers.'*'^ 

The  gas  temperature  as  determined  from  these  equations  specifies 
the  collisional  relaxation  rates  as  shown  in  Appendix  B.  In  addition  to 
these  rates  the  gas  temperature  determines  the  gas  flow  velocity  and  number 
densities  from  the  following  equations: 


v = (T/T  )(v) 
o o 


(D:  27) 


n = (T  IT) (n  ) 

a o xo 


(D: 28) 


where  the  o subscripts  denote  the  value  at  the  entrance  to  the  laser 
cavity,  i.e.,  at  x = 0.  These  expressions  indicate  an  important  difference 
between  constant  pressure  and  constant  volume  solutions  of  the  basic 
rate  Eqs . (3:1-3:71.  In  the  constant  volume  case  the  gas  flow  and  number 
densities  are  constants  whereas  in  the  constant  pressure  case  the  basic 
rate  equations  must  be  updated  by  Eqs.  (D:  27-28)  after  each  successful 
numerical  integration.  The  COFFEE  laser  calculations  assume  that  the 
pressure  is  constant. 
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APPENDIX  E 
STIMULATED  EMISSION 

Tne  interaction  between  10.6  ym  radiation  and  excited  CO 
molecules  leads  Co  stimulated  emission  and  absorption  between  the  00°1 
and  10  0 levels.  Stimulated  transitions  between  higher  lying  levels  are 
also  permitted  in  the  perfect  simple  harmonic  oscillator  approximation. 
However,  these  higher  levels  are  sparsely  populated  at  the  temperatures 
ol  interest,  and  the  enharmonic! ties  actually  present  in  the  CO  molecule 
Shift  the  wavelengths  of  these  transitions  into  regions  of  low  overlap 
between  the  linewidths.  Thus  only  transitions  between  the  00°1  and  10°0 
levels  need  be  considered.  Population  changes  due  to  photon-stimulated 
transitions  between  two  vibrational-rotational  levels  are  given  by  the 
relationships 


V ' Vll  - "yj'S2l' 


( E:  29) 


V ‘ 'Vi'  + \r< 

whire  S12  and  are  the  stimulated  absorption  and  emission  transition 

rates  respectively,  J and  J'  are  rotational  quantum  numbers  of  the  lower 
and  upper  states  respectively,  and  y denotes  either  the  P or  R vibrational- 
rotational  branch.  In  the  case  of  C02  the  zero  spins  of  the  identical 
nuclei  exclude  all  odd  values  of  J and  all  even  values  of  J'.  Thus 


J'  = J - 1 for 


Y = p 


J'  = J + 1 f 


(£30) 


or  y = R 
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‘T 


i. 


r r 


are  the  only  allowed  rotational  transitions  between  the  00°1  and  10°0 
vibrational  levels.  Due  to  the  changes  in  population  required  by 
1-q.  ( E:  29)  the  ASM  will  lose  k'9^  joules  and  the  SSM  will  gain  k.'0 
joules  of  energy  for  each  photon  emitted. 

ihe  population  densities  of  the  upper  and  lower  laser  levels 
can  be  written  in  terms  of  the  Boltzmann  factors  as 


nyj • = Ejt  a(l-a) (1-s) (1-b)  N£ 


V = s<1-a)U-s)(l-b)^  Nc 


( E31) 


where  is  the  fraction  of  molecules  in  the  Jth  rotational  level. 

This  fraction  ® is  given  by 

= 2(2J+l)(0r/T)e'J(J+1)0r/T  (a32) 

whenever  the  gas  temperature  is  much  greater  than  the  characteristic 
rotational  temperature  9f,  which  is  certainly  valid  for  C02.  The 
factor  of  2 is  the  nuclear  spin  symmetry  number.  Although  the 
characteristic  rotational  temperature  for  the  upper  state  is  lower  than 
the  value  for  the  lower  state,  the  difference  is  only  0.8%  and  can 
be  neglected.  ^ 

The  stimulated  emission  and  absorption  transition  rates  are 
related  by  the  principle  of  detailed  balance: 


_ _ 

QPJ 

s = 

r 2j~i 

sPJ 

12 

!_2  J+l  j 

S21 

. 

= 

~ 2 J-f-3  j 

qRJ 

b12 

_ 2J4-1  J 

S21 

(E: 33) 


where  Eqs.  ( E: 30)  have  also  baen  used.  By  combining  Eq.  (E:29)  through 
(E:33)  and  arranging  the  terms  in  the  proper  format,  the  stimulated 
transitions  in  the  basic  rate  equations  are  ob tained  where 
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r 


<.1  J 
yj  21 


PJ 


-2J  0 /T 

a - e 1 s / Q 


PJ 


RJ 


- (a-e2Utl,V\) 


RJ 


(E: 34) 


( E:  35) 


Qpj  = <2J-1) (20r/T)  e"J(J-1)Or/T 

( E: 36) 

Qrj  = (2J+3)(20r/T)  e-(J+2)U+l)&r/T 

These  equations  show  the  difference  between  the  P and  R rotational 
branches . 

The  stimulated  emission  rates  are  given  by  the  product  of  the 

photon  flu::  f^  at  the  appropriate  wavelength  and  the  stimulated  cross 

section,  o : 

yJ 


sn  ‘ fYJ  °tj 

(E: 37) 

lyJ  ' Vh'«o 

(E:  38) 

°yJ  ‘ 4vJ)A21/3rC^2 

(E: 39) 

W/nfe  *yJ  is  the  total  lj-8ht  intensity  at  the  laser  wave  number  v , 

(0)  o 

V is  the  llne  shaPe  factor  at  line  center,  and  is  the  spontaneous 

transition  rate  between  the  upper  and  lower  laser  levels.  The  spontaneous 

transition  rate  and  laser  wave  number  are  very  slightly  dependent50 

on  the  value  of  yj,  and  are  .201  sec"1  and  944.287  cm'1  respectively 

for  the  common  P(20)  transition.5^  5^ 
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At  low  pressures  the  line  shape  is  inhomogeneous  due  to 
Doppler  broadening  and  the  normalized  line  shape  function  is  Gaussian. 
The  value  at  line  center  is  given  by  the  relationships 


= A Jin  2 / 7t  B ^ 
u G 


$ : 40) 


B_  = v /8kT(£n2)/MWc 

b O 


(E:A1) 


where  is  the  Gaussian  full  linewidth  at  half  maximum  in  wave  numbers 
and  M is  the  CO  ^ mass  in  amu.  At  higher  pressures  the  line  is  homogeneously 
broadened  by  collisions  in  addition  to  the  homogeneous  natural  broadening 
of  the  transition.  In  this  case  the  value  of  the  line  shape  factor  at 
line  centor  is 


*<0)  - <2 /„>  >;} 


(E: 42) 


where  B is  the  Lorentzian  full  linewidth  at  half  maximum: 


■■  1 . 1_  v „ M /8kT/W 

*cttJ  i-  “V1  lV"“cO  i 

c,n,h,w 


(E:43) 


The  reduced  masses  p . are  easily  calculated,  and  the  optical 

cu2’1  52 
broadening  cross  sections  a have  been  measured  as  130,  91 

o 2 bu  2 * 1 

and  25  A for  collisions  with  C02>  N2  and  He  respectively.  The  optical 

broadening  cross  section  due  to  HO  has  not  been  measured,  but  can  be 

^ °2 

estimated  as  being  approximately  the  same  as  the  73  A measured  for  Ne, 

53 

which  has  a siailar  molecular  weight  The  natural  lifetime  in  Eq. 

( E: A3)  is  determined  by  the  sum  of  all  transition  rates  which  involve 
L 54 

the  upper  or  lower  laser  level: 
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T Ati  + S 2*?  + + R 

21  21  12  an 


(l+2An\ 

I-  ) 


+ R 


aw 


/ I+2Fw\ 

l ) 


+ R 


sb 


(^) 


ab 


1+2  Cb 
(l~b) 


31 


( E; 44) 


This  expression  does  not  include  electrical  excitation  which  is  usually 
negligible  relative  to  the  sum  of  the  other  transition  rates. 

Although  it  is  well  known  that  either  Gaussian  or  Lorentzian 
line  shapes  apply  in  most  practical  cases,  it  is  desirable  in  computational 
models  to  have  an  expression,  however  approximate,  which  can  be  used 
throughout  the  entire  pressure  range.  The  Voigt  profile55  accurately 
describes  the  combination  of  Gaussian  and  Lorentzian  profiles  but  is  too 
cumbersome  for  the  present  purposes.  Instead,  one  may  employ  a simple 
approximation  which  preserves  normalization  conditions  and  matches  both 
Gaussian  and  Lorentzian  expressions  at  the  appropriate  limits: 


'yd 


- U 


B 

C G 


+ *L  BL  )/(BG  + \ > 


which  reduces  at  line  center  to 


.(0) 

yJ 


/0)G  + <,> 


(E : 45) 


(E : 46) 


The  expression  can  be  used  at  all  pressures  and  temperature  to  yield 
the  stimulated  emission  transition  rates,  and  to  provide  a good 
estimate  of  the  Voigt  profile  over  the  transition  range  between  pure 
Gaussian  and  Lorentzian  line  shapes. 
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APPENDIX  F 


RADIATION  INTENSITY 

The  total  radiation  Intensity  I which  triggers  stimulated 

emission  ls  the  sum  of  the  spontaneous  radiation  intensity  lu  „ot 

collimated  by  the  cavity  mirrors,  and  the  spontaneous  plus  simulated 

radiation  intensity  i;,  which  is  collimated  by  the  cavity  mirrors. 
That  is, 


= TyJ  + XyJ 


(F:47) 


The  uncollimated  light  intensity  is  given  by 


I . - KP 
YJ  yj 


K - V // 

U 


(P.-48 


where  K is  the  power  density-to-light  intensity  conversion  factor  and 
YJ  15  thG  sP°ntaneo^  emission  power  density  within  the  gas: 


PyJ  = h'cv0  A2iQYJ(1~s)(1~b)2  a(l-a)N 


(F:4? 


The  geometric  parameters  used  in  constants  such  as  K are  easily  derlvab 
from  the  resonator  cavity  design,  and  are  given  in  Appendix  A. 

The  collimated  light  intensity  is  determined  by  Eq.  (3:7)  „hi< 
is  basically  the  equation  for  defining  gain.  The  last  term  is  the 
contribution  from  a constant  input  radiation  intensity  I to  be  amplify 
the  laser.  This  term  is  tero  if  the  laser  is  used  as  an  oscillator 

rather  than  an  amplifier.  The  light  inte„sity-to-p„„er  density  conversi 
factor  is  given  by 


K"  - A /V 
o c 


where  the  geometric 


parameters  are  given  in  Appendix  A. 


( F : 50) 
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Of  the  light  generated  by  spontaneous  emission  only  a very 
small  fraction  possesses  the  proper  spatial  mode  for  amplification  by 
the  resonator  cavity.  This  fraction  .is  given  by 


e , = >}|(?)  / (8ttv  2V  ) 
yJ  yJ  o c 


( F: 51) 


ihus,  the  second  term  on  the  right  side  of  Eq.  (3:7)  is  the  spontaneous 
emission  which  triggers  the  coherent  stimulated  emission  amplified  by 
the  laser  cavity. 

The  first  term  on  the  right-hand  side  of  Eq.  (3:7)  is  the  net 
gain  of  the  medium,  given  by  the  difference  between  the  total  gain  a 


and  the  average  cavity  loss  per  unit  length, 


YJ 


The  latter  term  is 


£,J  ' b,J  - ' nI (1'alYJ) <1'a2YJ) <1_TlvJ) a'T 2yJ) 1 /2L 

(F:52) 

where  b is  the  scattering  loss  per  unit  length  inside  the  gas,  L is 

the  distance  between  the  mirrors,  T is  the  transmissivity  and  a . 

iyJ  y iyJ 

is  the  total  loss  per  bounce  at  the  ith  mirror.  The  mirror  loss  term 
includes  scattering,  absorption,  and  diffraction.  The  total  gain  is 
obtained  from  the  product  of  the  stimulated  emission  cross  section  and 
population  inversion: 


a = a .An  , 
yJ  yJ  yJ 


The  population  inversion  is  determined  in  Appendix  E as 


(F:  53) 


AnyJ  “ VyJ  (1"s)(1_b)‘(1"a)  (*54) 

This  expression  is  positive  for  an  amplifying  medium  (although  the  gain 
may  not  be  above  threshold)  and  negative  for  an  absorbing  medium. 
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Although  the  net  gain  determines  the  radiation  Intensity 
within  the  laser  medium,  the  useful  output  Is  that  light  Intensity 
transmitted  by  the  end  mirrors,  namely 


IQut  = 2 Z (tiyj  + t2yJ)IyJ 

This  is  the  final  result  desired  from  any  laser 
is  necessary  because  at  each  mirror  only  half 
traveling  toward  the  mirror. 


(F:55) 

kinetic  model.  The  factor 
of  the  cavity  flux  is 
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